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ABSTRACT 


When a beam of extraordinarily polarized light 1s directed through a poled, single 
crystal of BaTiO, it fans in the direction of the optic axis. This thesis research 
investigates the origin of this fanning, or asymmetric self-defocusing, in barium 
titanate. Experimental observations and theoretical simulation suggest that this 
phenomenon is due to near forward stimulated photorefractive scattering. 
Additionally, observations of unique variegated beam patterns during phase 
conjugation are reported, including internally pumped rings and motion invariant 


patterns. 
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I. INTRODUCTION 


During the last decade the photorefractive crystal barium titanate (BaT10,) has 
received a large amount of attention within the scientific community. Barium titanate 
has long been of interest to acousticians and solid state physicists because of it's 
photoacoustic and piezoelectric properties, but only recently have optical physicists 
become very interested in this crystal. The cause of this sudden and intense interest 
was the discovery of photorefraction and optical phase conjugation. Originally found 
in liquids and gases [Ref. 1], optical phase conjugation has the ability to reconstruct a 
distorted wavefront (sometimes loosely referred to as time reversal). Barium titanate 1s 
one of many photorefractive crystals which can exhibit phase conjugation. 

There exist a number of intriguing optical properties of single crystals of barium 
titanate, but by far the most interesting 1s its ability to phase-conjugate at continuous 
wave low powers, with no external apparatus. This effect, known as self-pumped phase 
conjugation (SPPC), is the subject of extensive current investigation. The goal of this 
research is to determine the origin of self pumped phase conjugation in BaTiO}. 

Currently there exist two explanations of the cause of SPPC in BaTiO,. One 
theory, developed by Feinberg and his students at UCLA [Ref. 2], attributes the 
production of the phase conjugate reflection to the photorefractive effect via 
Degenerative Four Wave Mixing (DFWM). According to this theory, an asymmetric 
change in the index of refraction across the path of an incident beam causes the beam 
to bend. This bending eventually leads to a geometry which, through two internal 
reflections (often termed a CAT corner), allows the beam to intersect itself, leading to 
DFWM. The second theory, originating with Lam and his co-workers at Hughes 
Research Laboratory [Ref. 3], attributes SPPC to Simulated Photorefractive Scattering 
(SPS), or Two Wave Mixing, similar to that exhibited in gases and liquids with two 
frequency disparate input beams of very high powers. With the goal of proving one or 
the other of these theories, this thesis begins with a general background and theory of 
phase conjugation as it applies to BaTiO, (Chapter II). Chapter III describes PSIM 
(Photorefractive SIMulation), á simulation of the phenomenon of beam fanning due to 
a change in index of refraction within the incident beam path arising from the 


photorefractive effect. This simulation enables the prediction of the change in index of 


refraction and the index profile necessary to account for SPPC with only the 
photorefractive effect taken into account. 

Chapter IV introduces experiments to determine the effect of incident beam size 
on various crystals of BaTiO,. Chapter V summarizes the results and analyzes the 
data presented. The conclusion is reached that Stimulated Photorefractive Scattering is 
responsible for the process of beam fanning which leads to self pumped phase 
conjugation in BaTiO,. The actual origin of the phase conjugate beam is attributed to 
both SPS and DFWM. 

Additionally, Appendix A contains a summary of the work I did at Los Alamos 
National Laboratory under the direction of Dr. R. A. Fisher in February and March, 
1986. This work describes and catalogs internal beam production in barium titanate 
during SPPC. Some internal beam patterns not previously reported are introduced and 
explanations for these phenomena are proposed. The results of the experimental work 
in the body of this thesis have been submitted for publication to Optics Letters. An 
expanded version of Appendix A was presented at the XIV International Quantum 
Electronics Conference (1986), and has been submitted for publication in the Journal of 
the Optical Society of America B in collaboration with R. A Fisher and A. V. Nowak. 
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II. BACKGROUND AND THEORY 


A. INTRODUCTION 

Optical phase conjugation, first discovered in 1972 [Ref. 1], relies on a specific 
nonlinear property of a substance to create the formation of the complex conjugate of 
an incident wave. The reflection involved is not the familiar specular reflection seen in 
everyday life, but a reflection that retraces the path of the incident beam exactly, 
regardless of the direction or angle of incidence. The exact retracing of the incident 
beam path and the formation of the complex conjugate of the incident wavefront are 
the two unique and most important effects that make phase conjugation both 
interesting and useful. In barium titanate this phase conjugate beam 1s directly a result 


of the photorefractive effect. 


B. ORIGINS OF THE PHASE-CONJUGATE BEAM 

The origin of the phase conjugate beam in photorefractors is light which 1s Bragg 
scattered by an index grating formed by the light in the material. This process is 
known as wave mixing. There exists both two-wave mixing (scattering) and four-wave 
mixing in barium titanate. Both of these processes will be treated theoretically in 
general terms; the specifics of wave mixing in BaTiO, will be treated later. 

To understand the origins of phase conjugation one must begin with the wave 
equation for a homogeneous, isotropic, time-independent, nonmagnetic and 
nonconducting dielectric. The polarization term is divided into its linear (Pj) and 
nonlinear (P_,) parts, and the linear portion is absorbed into the linear permittivity of 


the medium, €. This leaves the wave equation as: 


3E  QP 
VE-zu—- ni 
BÉ H ĝt? 





(2.1) 


For ease of manipulation an approximation to this equation is made that reduces 
the second order nonlinear wave equation to a first order equation.! This 
approximation is derived by Fisher and Yariv [Ref. 5: pp. 9-11] and for it to be valid, 


the envelope of the pulse must not change appreciably during an optical period, hence 


Some calculations can be made without this approximation but they are quite 
tedious (e.g. [Ref. 4]). 
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the term Slowly Varying Envelope Approximation. In making the approximation it is 
necessary to assume that the electromagnetic waves under consideration are plane 
waves, and also to extract only the portion of the nonlinear polarization that 
propagates with the same frequency and wave vector as the modulation function of the 
electric field. This, referred to as the phase matched portion of the polarization, is 


given by 

P= P(z,t) exp+i(mt — kz), (2.2) 
in the wave given by 

E = E(z,t) exp+i(wt — kz). (2:3) 


Here the boldface letters (E and P) are the fields and the conventional letters (E and P) 
are the field envelopes, which vary little in an optical period and in an optical 
wavelength. All other terms are ignored because they are not phase matched and 
therefore cannot couple to the electric field. The Slowly Varying Envelope 
Approximation, or SVEA, is then invoked. That is: 


Ik2E] > > E > > ZE, (2.4) 
OZ oz? — 


Fisher and Yariv then show that the wave equation with SVEA becomes a first-order 


equation relating the envelope functions E and P: 
— + J/tu— = + P. (2.3) 


This resuit is the SVEA wave equation for a plane wave traveling in the +z direction. 
which shows that. given a nonlinear polarization, an electric field can be produced. 
This development can be generalized to the case of many waves present within a single 
medium, resulting in the origin of the nonlinear processes that can lead to phase 
conjugation. 


The polarization term can often be written as a power series expansion: 
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P = EX(E) = XWE + x0E? + x0ME3 +--> (2.6) 


where X(E) is the susceptibility of the material in which the wave is traveling. The 
expansion of the susceptibility in a power series separates the effects of the nonlinear 
material into distinct groups that categorize nonlinear effects [Ref. 5: p. 13]: 


e X(U. These are the linear properties that are the subject of classical optics and 
are not directly responsible for nonlinear effects. 


e XC). The second order effects are in general called three-wave mixing and 
include: second-harmonic generation, optical rectification, parametric mixing 
and the Pockels effect. 


e. XO). Third order effects include some of the most popular effects responsible 
for phase conjugation such as: third-harmonic generation, nondegenerate four- 
wave mixing, Raman scattering, ac Kerr effect (degenerate four-wave mixing), 
de Kerr effect, Brillouin scattering, two photon absorption, dc-induced 
harmonic generation and Stimulated Photorefractive Scattering. 


The two methods of creating a phase-conjugate reflection that are directly applicable to 
the experimentation described here are scattering and degenerate four wave mixing. 
l. Four-Wave Mixing 

In 1977 Hellwarth (Ref. 6] showed that it was possible to generate a time- 
reversed replica of any monochromatic-beam wave pattern. The method employed the 
interaction of an incident beam with counter-propagating pump waves in a 
homogeneous, transparent, nonlinear medium (see Figure 2.1). This method, called 
four-wave mixing, has become one of the most popular methods for phase conjugation 
because of the almost universal availability of the necessary equipment and proper 
media. 

The description of wave mixing begins with the nonlinear polarization derived 
from the energy function. As shown by Pepper and Yariv [Ref. 5: pp. 26-33], isolating 
only the third order term gives 


e. 4x3 
pr 4X l2pcF1E2E9 I (27) 

I consider here degenerate four-wave mixing (DFWM), where 
0,=0,=0,=0, represents frequencies of the four beams. The subscripts | and 2 
refer to two strong, precisely counterpropagating pump beams, the subscript p refers to 
the weak probe beam, and the subscript c refers to the phase-conjugate output wave. 


It is assumed that the following conditions are met: 


I5 


PHASE CONJUGATE SIGNAL BEAM 





Figure 2.1 Schematic of the process of four-wave mixing [Ref. 5: p.421]. 
a) intense (nondepleted) counterpropagating planewave pump beams of cqual 
intensity (k,, k,), 
b) a weak probe beam (K,), 


c) both the intense pump beams and the weak probe bcam are at the same 
frequency ©). 


Since the two pump waves are very strong compared to the probe beam they form a 


standing Wave. To obtain the steady-state solutions, we can set 
— =0. 20) 


So, assumung the SVEA and considering only the phase matched portion of the phase 


conjugate wave, then from Equation 2.5 


dE O0 JH 
x = dc P exp(ik.¢ r). (2.9) 
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The nondestructive buildup of the electric field amplitude can occur only if the product 
of the exponential phase factors in Equation 2.9 and P, is small (on the order of zero). 
If this is not so, then all of the successive contributions will not be in phase to add 


constructively. This total phase factor 1s: 
expl -i(k, +K)-k,-k,) e r]. (2.10) 


Since waves | and 2 are counterpropagating then, 
k, + k, = 0. (Con I) 


Although in practice this may not be exactly so because of nonlinear phase shifts due 
to unequal pump-wave intensities or non-aligned beams, it is very close [Ref. 5: p. 30]. 


Therefore the phase factor can be Zero only if 
k = —k . (2.12) 


Here it is evident that the nonlinearly generated field propagates in the opposite 
direction of the probe wave. 

Assuming that the factor in the exponential of Equation 2.10 is equal to zero, 
then combining Equation 2.7 and Equation 2.9 vields (within a multiplicative constant) 
(Ref. 5: p.36] 


OE. - 
io yn — X? EE ¿EN (2.13) 
Oz. " Ve 
This equation was solved concurrently by Yariv and Pepper [Ref. 7] and 
Bloom and Bjorklund [Ref. 8], and the solution was shown to be: 


2. dk] sin(IKIz) x cos(|k|(z-L)] 
E (2) -] TECOS cos(IKIL) Ee) 5: cos(IKIL) E (0) (2.14) 
cos(|k|z) _ « Sin{{k|(z-L)] _ x 
E — ———— E ——————— 
2) cos([k|L) dust Ik|cos([k|L) Ep > G 


Where 


* NI 


k" = @ X'BE, | (2.16) 


E (L) is the amplitude of the conjugate wave at its origin (the phase conjugate mirror), 
and E pL0) is the amplitude of the probe wave at its origin. 

In practice, the amplitude of the phase conjugate wave at its origin is zero (i.e. 
E (L)7 0), so that at the input point (z=0) of the probe beam 


x 


k i 
E (0) = iy enki END) (2.17) 


Equation 2.17 shows that at z=0 (the origin of the probe beam) the reflected field is 
proportional to the complex conjugate of the incident field. Therefore there exists a 
phase conjugate reflection. 
2. Scattering 
Phase conjugation by stimulated scattering was first reported by Zel'dovich, 
Popovichev, Ragulski and Faizullov in 1972 [Ref. 1]. As stated by Hellwarth 
[Ren ES 


In essence, the nonlinear polarization density, mediated by some driven (Raman-, 
Brillouin-, etc.) active excitation of the optical medium, couples normal-mode 
solutions of the linear Maxwell equation so as to create a new set of normal 
modes in backscattering with, in effect, complex propagation constants 
representing index change and gain. 


For example Stimulated Brillouin Scattering, or SBS, is a direct result of the 
electrostrictive effect seen in some materials. This effect is seen when the polarizability 
of a material is a function of pressure. In these materials it is possible for light to be 
scattered ov a sound wave allowing the coupling of a pair of light waves to a pressure 
wave. Another type oí scattering sometimes responsible for phase conjugation is 
Stimulated Raman Scattering or SRS. SRS is the scattering of light from molecular 
vibrations or rotations in a gas or liquid, or longitudinal-optical phonons in a solid. In 
SRS light waves are coupled because of polarizability changes with a molecular 


coordinate. Of particular importance to this research is Stimulated Photorefractive 
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Scattering in barium titanate. In SPS it is the photorefractive effect which 1s 
responsible for the production of a periodic index of refraction grating; this grating 
then Bragg scatters light in the backward direction which comprises the phase 
conjugate beam. 

The mathematics involved in all three types of scattering 1s similar and can be 
treated together. To this end I outline the work of Hellwarth [Ref. 5: pp. 177-180] and 
discuss the existence of the phase conjugate wave assuming scalar or symmetric 
scattering with a monochromatic, multimode, incident wave in a wave guide. The 
results are largely independent of the nature of the wave guide provided that a specific 
TEM mode that propagates in the forward direction also propagates in the backward 
direction. 


A wave E,(r) may be expressed as: 


A 


E,(r) = E A,e,, expli(k z — Dz)/2], (2.18) 


where f is the attenuation coefficient, € m 1$ the normalized transverse-mode pattern, | 
and the propagation constants may vary with wave angular frequency, that is k ne 
K AY): 

The electric field of the backscattered wave may then be written as 
Ref. 5: p.180] 


E, = E B e exp[-ik zfz/2 — yz/2), (2.19) 


where k, = k (©), and y is used to balance any nonlinear terms from the fields 
interaction with the medium. If y=0 in the above equation, then this is the solution 
to the linear Maxwell equations for a wave in a wave guide. 

If in Equation 2.19 y=0, then there is one configuration that allows the y of 
one of these modes to have a positive real part significantly greater than that of any 
other. This means that a single wave (reflection) may dominate the backscattered 


waves if the necessary conditions are met. The desired result is that 


E,(t) = nE, (r), (2.20) 
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where n is a constant. Substituting Equations 2.18 and 2.19 into Equation 2.20 gives 
the polarization density necessary to balance the extra terms in Maxwell’s equations 
that come from the last term in Equation 2.19 . According to Hellwarth [Ref. 9] the 


necessary conditions are satisfied if 


2 


(0 A FEA A 
Yk,(0)B, 7 n—7 T JJdxdy e, *ee, e KrniinA;A;B, , (2.21) 


where 
Kmin * jur) exp[AkLz]dz, (2.22) 
and 
Ak — k. (c) k(v)-k(V)-k, (o). (2.23) 


In Equation 2.21 I have assumed that the scattering considered is the pure scalar or 
symmetric type, reducing the fourth rank susceptibility tensor to a scalar quantity. 
This quantity, and all other constants, I have absorbed into n. 

The necessary conditions may be determined for special cases, either by exact 
solution or perturbation methods developed by Hellwarth. These solutions indicate 
that a backscattered solution that is nearly phase conjugate predominates because it 
has gain nearly twice that of any other mode, provided that [Ref. 5: p. 177]: 
the interaction length is not too long, 


the Stokes shift (v — @) is not too large, 


og 


the number of guide modes falls within a certain large range, and 
d. the total number of guided modes is not too large. 

The solutions that led to these conditions. presented by Hellwarth in 1978 [Ref. 91, 
unfortunately do not take into account competing nonlinear effects, and thereforemme 
good only to the extent that they predict the major effects observed in experimentation. 
It is observed that the backscattered wave is not always an exact phase conjugate of 
the input wave, possibly due to some of these competing effects. These analyses apply 
equally as well to SPS, SBS and SRS [Ref. 3]. 


C. DISTORTION CORRECTION BY PHASE CONJUGATION 

Although not immediately applicable to this investigation, the ability of a phase 
conjugator to correct for phase distortions is one of the primary aspects driving 
research in this area. As a practical point, this ability is an easy way to determine that 
any given observed beam is a phase conjugate of another. 

It was shown above that, at least for the two cases considered, a wave can be 
generated propagating directionally opposite to an input wave, with it’s envelope 
function proportional to the complex conjugate of the input beam’s envelope function. 
This phase conjugate wave may be considered a separate wave that must satisfy the 
same wave equation as the input wave. This means that if the wavefront of the input 
wave is distorted in some way, then the phase conjugate wave will propagate 
backwards along the same path and its envelope function will remain everywhere the 
complex conjugate of the input wave's; that is, its wave fronts will coincide with the 
input wave fronts at every point. Thus a distorted wave front will be restored to its 
original shape as the phase conjugate wave traverses the distorting medium in the 
reverse direction. The proof is quite simple, as shown by Fisher and Yariv [Ref. 5: p. 
17], and is accomplished by complex conjugating the wave equation for the forward 
going wave. The result is the equation for the conjugate wave propagating in the 
opposite direction. So a wave traveling in the opposite direction of the incident wave, 
with an envelope function the complex conjugate of such a wave, satisfies the same 
wave equation as the incident wave. This means that the backward going wave will 
remain everywhere the complex conjugate of the incident wave, and therefore a 
distorted plane wave will be restored to its original configuration when it is passed back 
along the same distorting path. This distortion correction capability is the most visible 
and fascinating effect attributable to phase conjugation at the present time. 

An incomplete list of some applications of phase conjugation with emphasis on 
those currently using barium titanate as the phase-conjugator is [Ref. 5,10,11]: 


a. Holography 


En 


Distortion correction 

Laser oscillators 

Interferometry 

Laser beamsteering 

Associative memory investigations 


Image amplification 


S MEAL 


Laser coupling 
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i. Image addition, subtraction and enhancement 
j. Optical bistability 
k. Automatic pointing and tracking. 


D. PHOTOREFRACTORS AND FOUR WAVE MIXING 
In photorefractive materials the index of refraction changes with incident optical 

energy. This effect was first discovered in the 1960's [Ref. 12,13] and has since led to 
the development of optical recording sensitivities comparable to that of the silver halide 
emulsions used in photography [Ref. 14]. The process involved in the effect is 
summarized by Feinberg [Ref. 5: p. 418] as follows: 

a. Light causes charge to migrate and separate in a crystalline material. 

b. The separation of charge produces a strong electrostatic field. 


c. The electrostatic field causes a change in the refractive index of the crystal by 
the linear electro-optic effect.’ 


In the photorefractive effect the important factor in the change in index of 
refraction is not the intensity of the incident light, but the total energy of the incident 
radiation; the power only determines the speed of the reaction. Also because of the 
type of effect concerned, the actual result is dependent on the relative intensity of the 
incident radiation. 

Consider a thin beam of light incident on a photorefractive crystal that is 
otherwise dark. The light will liberate charge carriers from the lighted area of the 
crystal and these carriers will then settle back into a place within the crystal structure. 
If the carriers settle back in the area of the crystal that is lighted they will again be 
liberated. If, however, the carriers land in a dark area of the crystal, they will stay 
where they land due to the small dark conductivity of the crystal. Eventually there will 
be few charge carriers in the lighted area and many in the dark area. This 
disproportionate distribution creates an electric field between the dark area and the 
lighted area that creates a change in the index of refraction of the crystal through the 
linear electro-optic effect. 

This process of migration has been aptly described by RK. Al Fisher sass 
ilyswatter effect (Ref. 16]. It is as if the charge carriers were flies at a picnic table and 


the light was a flyswatter suspended from above, slowly swinging across the table top. 


The linear electro-optic effect is the change in the index of refraction caused by 
(and proportional to) an applied dc electric field. This effect 1s seen in crystals that do 
not exhibit inversion symmetry, such as barium titanate in the tetragonal phase 
[en 15 p 225]. 
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Each time a fly lands on the table within the area of the flyswatter’s reach, he will 
immediately fly up as the swatter returns toward him. If, however, the fly lands on the 
table in an area away from the swatter he will be left alone and will stay on the potato 
salad. In time therefore, there will be no flies on the table in the area where the 
flyswatter can reach and many where it cannot. This process effectively pumps the 
flies to regions on the table which are away from flyswatter. 

The model of charge transport described above is known as the hopping model 
and was first proposed by Chen to explain the change in refractive indices in LiNbO, 
and LiTaO, [Ref. 17. There also exists a diffusion model proposed by Amodei 
[Ref. 18] that may dominate when the spacing between light and dark areas within the 
crystal is less than one micrometer. 

From the flyswatter analogy it is clear that what is important is the optical 
intensity relative to the ambient. This entire explanation assumes the existence of 
charge carriers in a crystalline material. In most materials exhibiting the 
photorefractive effect the origin of the charges is unknown, but it is assumed that they 
inhabit low-lying traps formed by impurity sites in the crystal. For example, in 
BaliO,, Fe?* and Fe?* impurities may create the charge carriers as shown by Klein 
and Schwartz [Ref. 19]. It has also been proposed, and there is extensive evidence to 
support the theory, that oxygen vacancies are the donors and barium vacancies are the 
acceptors [Ref. 20]. Whatever the dominant species is, it is assumed that there are a 
large number of empty receptor sites in the material available to the liberated charge 
Carriers. 

If, instead of a single beam, two beams of identical wavelength and equal 
intensities (one called a pump wave and one called a probe wave) intersect within a 
photorefractive crystal, the interference pattern they form will cause a periodic optical 
field to be set up in the region of intersection, resulting in a periodic charge separation. 
Through the electro-optic effect, a periodic variation in index of refraction will result. 
This index grating is then able to Bragg scatter a pump wave in the opposite direction 
to produce a phase conjugate beam. The formation of the index grating is shown in 
Mente 2.2, adapted from Feinberg [Ref. 5: p. 422]. 

It is important to note that photorefractive phenomena are not merely four-wave 
mixing phenomena. Four-wave mixing is a consequence of an index change at some 
given point at some given time (x,t). Photorefractive phenomena however, are 


dependent on the history of the medium and are not spatially local effects. That is, the 
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Figure 2,2 Index grating formation due to the photorefractive effect (Ref. 5]. 


effects at (X,t) are a consequence of actions at (x-óx, t-ót). This is clearly seen in 
Figure 2.2: 

The top curve of Figure 2.2 shows the intensity as a function of the position 
along the path of the beam. The second curve is the charge density due to this 
intensity variation. The electric field intensity 1s presented in the third curve (note the 
important on period shift relative to the charge). Finally the variation of the actual 
index of refraction within the crystal is shown in the last curve. 

The spatial shift between the change in index of refraction and the variation in 
intensity due to the interference pattern (one fourth of the grating period) is an 
important effect. This happens because the peak of the charge distribution is a center 
of right-left symmetry, and theretore there can be no clectric tield at that point. This 
shift leads to the ability to transter power from one beam to another, a process known 
as two-beam coupling. Two-beam coupling, first reported by Staebler and Amodei 
(Ref. 21], is due to the scattering of one beam into another by Bragg diffraction as they 
pass through the grating they have mutually produced. The result is that one beam 


experiences gain at the expense of the other. 


Do 


One interesting effect of four-wave mixing in photorefractors 1s seen in the phase 
conjugate mirror with continuous wave gain. This was first demonstrated by Feinberg 


in 1980 [Ref. 22] using the configuration shown in Figure 2.3 . 


reflective 
surface 
BaTiO, 
crystal 


pumping beam pumping beam 


output beam 


Figure 2.3 Formation of an optical cavity with a phase conjugate mirror [Ref. 5]. 


In this configuration a phase conjugate mirror (PCM) and an ordinary nurror (Ml) 
can form an optical cavity in which a beam will oscillate. One of the more dramatic 
demonstrations of this type of oscillation is reported by Feinberg [Ref. 5: p. 436] in a 
configuration using an ordinary kitchen spatula as one end of the cavity. Just as in 
any conventional laser oscillator, the oscillations can begin with noise (as little as a 
single photon) that makes one successful round trip within the cavity. With each 
additional round trip the system sees gain due to two-beam coupling and in a sliort 
period of time a visible CW beam may be observed between the mirror and the phase- 


conjugator. 


E. STIMULATED PHOTOREFRACTIVE SCATTERING 

The second theory which accounts for phase conjugation in BaTiO, was 
Originally proposed by Lam in 1985 [Ref. 5]. Lam proposed that the origin of the 
phase conjugate beam is two wave mixing or scattering. This process of SPS begins 
with random inhomogeneities in the photorefractive medium which give rise to noise 


photons from the Rayleigh scattering of the input beam. In the words of Lam [Ref. 3]: 


The noise photons propagating along with [sic] axis of the input beam undergo 
parametric scattering, resulting in an unidirectional gain. The phase conjugate 
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wave arises from the component of the scattered wave which experiences the 
maximum gain. 


Later the same year Chang and Hellwarth conclusively demonstrated SPS in BaTiO, 
(Ref. 23]. To do this they had to index match the surfaces of the crystal to make 
accessable angles which would not produce the characteristic corner reflection (i.e. the 
CAT corner) producing DFWM. 

Later, Vallev [Ref. 24] defined SPS as the selective amplification through the 
photorefractive effect of optical radiation scattered by medium inhomogeneities. Thus, 
SPS is a two-wave mixing process in which the incident beam interferes with scattered 
light to create the necessary grating through the photorefractive effect. As soon as the 
gain for an individual noise beam exceeds its loss, then that beam will become 
dominant, increasing the grating strength and hence the backscattered beam gain. This 
process will continue until a steady state value is reached. 

The mathematical equations, as derived by Lam, are of the identical form to 
those encountered in stimulated scattering in a waveguide (presented above), so the 


analysis of phase conjugation by SPS has already been addressed. 


F. BARIUM TITANATE 
1. The Crystal 

The experiments reported here were performed with one of the more popular 
photorefractive crystals capable of creating a phase conjugate reflection. Above the 
Curie temperature of 120°C, BaTiO, is a cubic crystal and therefore centrosymmetrical 
and non-photorefractive. Between 5°C and 120°C however, the crystal is in the polar 
phase and has a tetragonal symmetry. The axis of fourfold rotation is the optical axis 
or c-axis. In the tetragonal phase barium titanate is ferroelectric, photorefractive, 
photoconducting, birefringent, electro-optic, pyroelectric, piezoelectric, photoacoustic 
and (slightly) electrostrictive. The manv possible effects make the crystal itself very 
difficult to understand, since often the results of the various effects cannot be 
separated, All results reported here involve BaTiO, tn the tetragonal phase. 

Since BaTiO, is chemically and mechanically stable at room temperature, it ts 
used extensively in many capacities and is therefore well studied. Jona and Shirane 
[Ref. 25: pp. 108-215] report extensively on the structure, properties and the effects of 
temperature, stress and electric fields on barium titanate in the tetragonal phase. 


Although at the time of printing the photorefractive capabilities were unknown, their 
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data represent a comprehensive picture of the crystal and forms the basis for analysis 
of the solid state aspects of barium titanate. More recent but less comprehensive data, 
emphasizing the optical aspects of BaTiO,, have been compiled by Wemple, et. al. 
[Ref. 26]. 

When used as a photorefractor, it is necessary that a single crystal of BaTiO}. 
be used and that the crystal be poled into a single domain. Poling is accomplished by 
cooling the crystal through the Curie point while an electric field on the order of 10° 
V/cm is applied across the optical axis. This poling of the crystal results in an electro- 
optic tensor which has a very large [42] element. Since the photorefractive effect uses 
the linear electro-optic effect, it is important to maximize this [42] element. This 
maximization is accomplished by using extraordinarily polarized light incident at an 
angle between zero and 90° to the c-axis. 

2. Beam Fanning 

When a beam of extraordinarily polarized light is incident on a poled, single 
crystal of barium titanate, it is seen to asymmetrically defocus in the direction of the c- 
axis prior to the onset of phase conjugation. This process, called beam fanning, was 
first discovered by Feinberg in 1981 [Ref. 27]. Beam fanning is clearly visible in Figure 
2.4 (the optic-axis is designated by the arrow). This fanning can be explained by either 
an asymmetric change in index of refraction within the incident beam path, or by 
stimulated scattering. 

a. Deflection by Photorefraction 

Beam fanning was proposed to be due to an asymmetric change in the 
index of refraction caused by the photorefractive effect. This process is represented in 
the curves of Figure 2.5 (from Feinberg [Ref. 27] ). The top curve shows the intensity 
curve of a Gaussian beam. The second curve shows the net charge due to 
displacement of charge carriers by the incident light. The third curve shows the 
induced electric field, and the last curve shows the asymmetric change in the index of 
refraction. It is assumed that at the center, and most intense portion of the beam, the 
gradient of the index of refraction is greatest. Therefore the center portion of the beam 
creates the observed effect, and the fanning is induced beginning at that point. In 
other words, beam fanning begins in the center of the beam. This is one of the major 


discriminating factors and will be referred to later. The second major point to consider 


?Note that even if the intensity profile of the incident beam is not smooth, the 


top curve of Figure 2.5 represents the envelope that the beam intensity must fall 
within. 
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Figure 2.4 BaTiO, exhibiting SPPC. 


about this theory is that as the diameter of the beam decreases the electric field 
increases, since the linear distance between charge distributions decreases. Since the 
photorefractive effect is proportional to the electric field, then beam fanning must 
increase with a decrease in beam diameter. This is clearly seen by examining the 
equations describing the electric field and the subsequent refractive index change. 

The electrostatic field induced by an incident beam in a photorefractor 1s 
given by [Ref. 27]: 


-2kgT(Xx + yy) 


2 
eNO, 


E(x) = 


expl(x* + y) 0,4, (2.24) 


where e is the charge or the charge carrier, kp is Boltzmann’s constant, I is the 
temperature of the crystal, and «,, is the incident beam diameter. The change in index 
of refraction in a crystal of BaTiO, for an extraordinary ray due to this field is 


approximately given by 
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Figure 2.5 The process of beam fanning due to photorefraction 
for a single incident beam [Ref. 27]. 


n(y) = -[E(x)*y] m ovn Tano), 5 


cosÜ 
n 
where n(0) = (me cosĝ + nes singy 1/2 , 9 is the angle the incident beam makes with 
the non c-axis of the crystal, Du is the index of refraction for an ordinary ray, D 1s the 
index of refraction of an extraordinary ray, and all terms that do not include the r}, 
element of the electro-optic tensor have bcen dropped, since its value exceeds all other 
elements by an order of magnitude. 

Inspection of these equations clearly shows that the change in index of 
refraction is inversely proportional to the diameter of the incident beam. 

b. Near Forward-SPS 

In 1985 Lam proposed that beam fanning is due to an effective nonlinear 
index of refraction which causes the crystal to behave like a diverging Iens [Ref. 3]. 
This nonlinear index of refraction (~ 107) is due to an intensity dependent frequency 


shift of the incident beam due to the bulk photovoltaic effect. Recently Valley 


D 


[Ref. 24] has suggested that beam fanning is actually near forward scattering due to 
SPS, in an extension of Lam's theory. 

As in other types of scattering, SPS is the selective amplification of 
radiation which is scattered by inhomogeneities in the medium. The amplification is 


accomplished through the photorefractive effect. This is schematically diagrammed in 
Figure 2.6. 





SCATTERED 
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Figure 2.6 The process or Sis | Rete.) 


In a crystal where one charge carrier dominates, the gain-length product 


per absorbed energy per volume at turn on is given by [Ref. 24] 


= At opr enteE slope (2:36) 
s Es E] 


t 


where r,rp is the effective Pockels coefficient for extraordinarily polarized light [Ret. 25], 
A = 2t/X, e is the carrier charge, t is the space charge field and refractive index grating 


decay rate, € is the dielectric constant (which is anisotropic), and 


EE, = kpTA ute, (2.27) 
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where kg is Boltzmann's constant, H is the carrier mobility (which is also anisotropic) 
and T is the temperature of the crystal. Finally, lp is the effective interaction length 


which Valley models with 
leg 7 L(ül-exp[-G,/(LsinO)p, —. (2.28) 


where L is the length of the crystal in the direction bisecting the scattered and incident 


beam, 0, is the incident beam diameter, and 


9 = (0, - 0,)/2. (2.29) 


The important item to note is that the gain-length product is proportional 
to the effective interaction length given by Equation 2.28 . Therefore the amount of 
Scattering due to SPS (both forward and backward) is exponentially dependent on the 
incident beam diameter in such a way that an increase in beam diameter will increase 
the scattering, and a decrease in beam diameter will decrease the scattering; this 1s 
opposite from that expected by the theory discussed in Section F.2.a above. Also, 
since SPS is dependent on the interaction of random noise photons with the incident 
beam, one would expect that beam fanning (near forward-scattering) could commence 
at the periphery of the beam in addition to the center. This scattering of energy should 
be more obvious at the edges of the incident beam since the intensity is lower there. 
Also, any light scattered toward the c-axis from inside the beam would take power 
from the periphery through two beam coupling, thus adding to the loss of intensity 
toward the beam edge. Considering these points, if the intensity profile of the incident 
beam were to be monitored during the fanning process, one should expect to see the 
intensity diminish from the periphery toward the center. Again this is opposite the 
effect expected from the prior theory. 

3. Self Pumped Phase Conjugation 
Beam fanning is intricately invoived with seif-pumped phase conjugation 
(SPPC). SPPC in BaTiO, was iirst discovered by Feinberg soon after the discovery of 
beam fanning [Ref. 2]. According to Feinberg, the process of beam fanning (due to an 
asymmetric change in index of refraction within the beam path) may result in a portion 
of the incident beam being bent in such a way that it is internally reflected off the 


corner of the crystal (see Figure 2.4). This reflected beam then intersects the incident 
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beam causing a steady state interference pattern. The light and dark areas of the 
interference pattern cause an index grating (through the photorefractive effect) which 
then Bragg scatters the incident beam in the process of DFWM, creating a phase 
conjugate beam. Lam alternatively theorized that the origin of the phase conjugate 
beam is two-wave mixing (SPS) as described above, and that auxiliary beams were not 
essential to the process. 

Of interest also is the time for phase conjugation to begin in barium titanate. 
The time to produce self-pumped phase conjugation can vary from a few seconds to 
tens of minutes depending on the intensity of the incident beam and the configuration 
used. Because of the typically large time for the onset of phase conjugation in BaTiO,, 
there is extensive work presently underway to try to decrease it. Some work on 
changing crystals after growth has been reported [Ref. 20], but the present emphasis 1s 
on the doping of the crystals during growth [Ref. 29]. 

The exact processes involved in beam fanning and phase conjugation are still 
unresolved in the literature. In the following chapters I will examine these theories as 
they apply to observed effects and, since in the case of beam fanning there are 
mutually exclusive indications, I will show that beam fanning must be due solely to 


near forward SPS. 
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IM. PSIM: A SIMULATION OF THE PHOTOREFRACTIVE EFFECT 
IN BARIUM TITANATE 


A. INTRODUCTION 

To investigate the photorefractive effect in barium titanate, I designed a 
computer simulation which plotted the internal beam pattern within the crystal. This 
simulation, entitled PSIM (Photorefractive SiMulation) only accounts for the 
nonlinear photorefractive effect, not for SPS or other competing nonlinear effects. In 
so doing, I isolated the photorefractive effect from other effects and determined the 
necessary change in index of refraction across the beam path to account for beam 
fanning which is typically observed in the laboratory. The source code for PSIM 1s 
found in Appendix B. 


B. THEORY, 

As discussed above, the photorefractive effect may be solely responsible for the 
observance of beam fanning, which ts in turn responsible for self-pumped phase 
conjugation [Ref. 2]. The change in index of refraction of a given crystal, An, may be 
given by the equation [Ref. 14] 


An(z)= — (1/,)n?re¢p El2), (3.1) 


where n is the index of refraction, rypp is the effective Pockels coefficient and E(z) is the 
space charge electric field. 

Values for rfp range generally in the picometers per volt range. Of seven crystals 
used by Klein and Schwartz [Ref. 19], rep ranged from 4.2 x 10717 to. 12.0 x 107? 
meters per volt. Although these values are very small, the space charge electric field 
may be very high, on the order of 10° volts per meter [Ref. 5: p. 418]. Using these 
values it becomes obvious that the change in index of refraction due to the 
photorefractive effect in barium titanate may be as high as 107^ . Though not an exact 
number for any given crystal, this may be accepted as close to an upper bound for An 
(denoted here as An nax) for most of the crystals used. Günter [Ref. 11: p. 226] gives 
Bun as 2.2 X 107, limited by the possible space charge field. 
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From the above simple calculations it may be assumed that for the 
photorefractive effect to be solely responsible for self-pumped phase conjugation, the 
dimensions of any given crystal used must be such that a change in index of refraction 
on the order of 10°* is sufficient to cause the observed beam fanning. However, as 
previously mentioned, the gradient of the electric field and hence the index gradient is 
the important quantity relating to the photorefractive effect. Therefore, the beam size 
relative to the length of the crystal along the direction of propagation and the intensity 
profile must also play an important role in SPPC. These factors, in addition to the 
incident angle, were eventually incorporated into the simulation described here. The 
resulting simulation was used to determine the minimum change in index of refraction 


necessary to account for observed beam fanning in any given configuration. 


C. THE SIMULATION 
1. Program Logic 

The simulation is written in HP basic and designed to run on a Hewlett 
Packard 300 model computer with at least two megabytes of random access memory. 
This two dimensional model takes into account only the change in index of refraction 
across the incident beam path due to the photorefractive effect. The output is a 
graphic representation of the eventual steady state beam geometry. 

The model uses a brute force approach to calculate the observed effects which, 
although time consuming and memory intensive, uses a simple algorithm. The 
program uses onlv Snell's law and Huygen's principle to calculate the internal beam 
pattern. Input into the program is the range of the change of index of refraction 
across the input beam, An, the crystal and beam size (in arbitrary units), the point of 
entrance of the beam into the crystal, and the angle of incidence of the incident beam 
relative to the normal. The program computes the mean index of refraction between 
two input extremes which then corresponds to the normal index of refraction (n~ 2.4). 


The important factor in the program is the difference between the two extremes, 


Animax - Logically it is this entity which drives the photorefractive effect. Note that 
ihe simulation is not concerned with how An... is created. 


Once the necessary parameters are entered, the internal angle of the incident 
beam is calculated and each point on the wavefront is then assigned an index of 
refraction value that it can impart to the crystal at its location. These points are then 


grouped into pieces of the wavefront (one to five separate pieces are normally used). 
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The assigned index value is dependent on An pax and an assigned profile which may be 
varied (more will be said about this profile later). The indices of each point of the 
wave front are then propagated in a straight line through the crystal, at which time the 
value of the index of each point in the wavefront is imparted to the corresponding 
point within the crystal. After completing this process the wavefront is propagated 
along the same path, allowing each point to move at a speed proportional to the index 
at the point in the crystal in which it finds itself. After a short time the process of 
wavefront movement is stopped and the direction of each: piece of the wavefront is 
calculated using Huygen’s principle. The original indices of the points within any given 
piece are then propagated in that new direction, imparting those values along the new 
path of the piece of the wavefront and the process begins again. The result of many 
iterations of this algorithm is a graphic output of the steady state internal beam 
pattern due to the photorefractive effect. 
2. Limitations 

There are several limitations to the model that must be discussed prior to the 
introduction of results. The first, and major, limitation is that of necessity each point 
on the wavefront cannot be treated individually. The wavefront must be treated in 
sections to determine an appropriate direction for a section and this makes the model 
only an approximation to the actual continuous case. The choice of the number of 
sections is primarily dependent on the chosen index profile. A profile similar to that 
shown in Figure 3.1 may be treated as a single piece since the direction of each 
individual piece will be the same as the direction of the entire wavefront. By contrast 
compare the profile of Figure 3.2. This profile appears to lend itself to the three pieces 
shown; however, in using only three pieces the effect around the extremum are 
minimized and it may be wise to consider the wave in five pieces. 

The difference in considering three as opposed to five pieces in the wavefront 
may be seen in Figures 3.3 and 3.4. In Figure 3.3 the index profile shown in Figure 3.2 
was used in three pieces. Note the distinctness of the three discrete sections. Figure 
3.4 shows the result of the same profile divided into five sections. It is still easy to 
mote the five individual sections under consideration, however the output :s noticeably 
different. Ideally it would be desirable to have a very large number of pieces used in 
wavefront calculations. However, the larger the number of pieces, the larger the 
number of points in the wavefront, and consequently the larger the beam size. To keep 


the beam/crystal ratio in the proper perspective it is necessary to increase the size of 
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Figure 3.1 Linear index profile considered is one piece. 


the crystal whenever increasing the size of the beam, and the maximum crystal size 1s 
limited by the available memory. The machine used for this research was limited to 
two megabytes of random access memory, Waich allowed for a maximum crystal size of 
500 units. To keep appropriate beam/crystal ratios, the beam size used (generally) was 
between 20 and 40 units and the number of pieces in the wavefront ranged from one to 
eight, depending on the index profile. 

A second limiting factor in the use of this simulation is seen in cases where 
there is beam crossover. This happens when the index profile has a larger gradient on 
the side away from the c-axis than it does on the c-axis side. An example of this type 
of profile can be seen in Figure 3.3. In this case the bedin intensity ap mo 
considered to be Gaussian and the more intense section of ¿he Deam (rearme 
nearest the center) dominates over the weaker section. This ignores any consideration 
of two beam coupling which would actually occur. 

The final major limiting factor is the inability of the model to fully account for 
a piece of the wavefront which breaks away from the main poruon of the incident 


beam. A case where this happens is seen in Figure 3.3 . In this case the indexmes 
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Figure 3.2 Nonlinear profile with possible divisions. 


refraction gradient across the piece is the same as when it was part of the entire beam. 
This may be a proper assumption if the beam was continuous during fanning; however, 
as can be seen in Figure 2.4 and many others presented here, there are often portions 
of the beam which leave the major portion of the incident beam entirely. Under the 
condition of a complete disconnection it is possible that the total index gradient should 


be based on An 


wavefront piece was within the incident beam; this was not programmed iuto the 


max and not on the index gradient as it was calculated when the 


modei. There are presently no data available indicating the effect of a detached piece 
of a beam in contrast to the continuous case. This problem is circumvented when a 


profile such as that shown in Figure 3.1 is used so that only one wavefront piece is 


considered. 


D. RESULTS 


As noted above, the output of the simulation 1s dependent on the index profile, 
the beam/crystal ratio, the angle of incidence and the position on the crystal where the 


incident beam enters. The criteria for success (of possible SPPC) of a given 
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Figure 3.3 PSIM output for profile shown in Figure 3.2. 


configuration was taken to be the ability of a portion of the beam to reflect off the 
crystal corner in such a way as to reintersect itself. This would allow for. the 
production of auxiliary beams which are a prerequisite to SPPC. Thercfore, a result 
such as that seen in Figure 3.5 would be considered capable of SPPC, and one such as 
that seen in Figure 3.6 would not. 

Since the purpose of this investigation is to model beam fanning, the most 
important parameter is the index profile. To concentrate on this parameter, a 
oeam/crystal ratio and angle of incidence that matched one of the configurations seen 
in the laboratory were chosen. The configuration chosen was an often used 
configuration with a beam/crystal ratio of 0.064 and an incident angle of 40? . The 
incident beam entered the crystal at a point E way across the face of the crystal. This 
corresponds to a crystal observed in the laboratory of 5.0 mm on a side with a beam of 


size 0.32 mm, incident 3 mm from the front side of the crystal. 
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Figure 3.4 PSIM output for profile shown in Figure 3.2 ; 
but broken into five separate pieces. 


` i 
` . 


Modeling was run with three index profiles. The rationale behind the use of the 
profiles will be explained in Section E of this chapter. In each case the values of 
Anmax Which caused successful SPPC due to beam fanning are presented. 

Following the theory of Feinberg (Ref. 2], the first index profile used was that 
shown in Figure 2.5. Successful SPPC from beam fanning was seen (Figure 3.8), 
however, note the beam fanning in the direction opposite the c-axis: this is never 
Bissesved in practice. 

The second beam profile used 1s shown in Figure 3.1. A typical successful run is 
shown in Figure 3.9. Although this approximates some observations of beam fanning, 


note the lack of a continued incident beam as one often (but not always) observes. 
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Figure 3.5 Beam fanning resulting in successful SPPC. GN 


The final profile used is shown in Figure 3.7. SPPC due to beam fanning was 
exhibited as is seen in Figure 3.10. It was this configuration which exhibited results 


closest to those observed in the laboratory. 


E. DISCUSSION 

The first index profile is taken directly from Feinberg [Ref. 2]. To arrive at a 
profile such as this it is necessary to assume that charge migration in the 
photoreíractor ıs symmetric, resulting in a net charge distribution such as shown in 
Figure 2.5. The beam pattern predicted by the simulation is shown in Figure 3.8. The 
major problem with this profile is the fanning of the beam in the direction opposite the 
c-axis due to the tail of the index profile. In his original work Feinberg only considers 
the central region since it is the most intense part of the beam, ignoring the tails of the 


profile. As can be seen from the model, the tails of the profile cannot be totaliv 
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Figure 3.6 Beam fanning not resulting in successful SPPC. 


ignored. To arrive at the observed effect with the index profile of Figure 2.5 it ts 
necessary to assume that barium titanate exhibits asymmetric self-focusing (i.e. self- 
focusing only in the direction of the c-axis). If asymmetric self-focusing was present, 
then the (less intense) light subject to the profile of the tails would be unable to fan 
against the c-axis. To date there has been no reported evidence of self-focusing of any 
kind in DaTiO,; however, this does not preclude the possibility since it would only 
affect a weak portion of the beam and may not be obvious to the casual observer. 
Note that the fanning against the axis on the left-hand side of the beam would be 
abated by two beam coupling and might not be readily discerned. 

An alternative explanation is that the crystal response to the beam is not 
characterized by the curves shown in Figure 2.5. Instead, all available charge may be 
totally excluded from all portions of the beam, resulting in the profiles shown in Figure 
STI. 


By 
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Figure 3.7  [ndex profile for simulation of beam fanning in BaTiO.. 
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In an effort to approximate an index profile which would exhibit the necessary 
behavior, the profile shown in Figure 3.1 was used. This beam profile was arrived at 
by making the important assumption that there is a total exclusion of charge carriers 
within the illuminated section of the crystal. This is an easy profile to work with since 
the outcome is independent of the nuniber of pieces considered in the wavefront, as 
previously mentioned. The result of the simulation with this profile is shown in Figure 
3.9. However, there are two major problems with this profile. The first one 1s the lack 
of extension of the incident beam as mentioned above. Note that in most of the 
photographs presented in Appendix A there is an extension of the incident beam. This 
is almost universally, though not exclusively, true (see Figure A.5 for example). 
Second, it is unrealistic in its sharp transitions. Both of these problems are solved with 
the introduction of the third index profile. 

The final index profile (Figure 3.7) resulted in a beam geometry which exhibited 
all of the requisite properties. This profile was arrived at by the rounding of the 
discontinuities of the prior profile; specifically, a sine function was used to arrive at the 


configuration shown in Figure 3.7. The results of the simulation are presented in 
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Figure 3.8 PSIM output using the index profile of Figure 2.5 


Figure 3.10. To arrive at this profile 1t would be necessary for diffusion of the charge 
carriers almost completely out of the incident beam path. The eventual diffusion of 
charge carriers completely out of the path may actually be expected since any light of 
intensity above ambient will cause any charge carriers to eventually be totally 
liberated.* The graphs of Figure 3.11 show the necessary charge, electric field and 


intensity profile of the incident beam necessary to arrive at an index profile similar to 
[t is assumed that the continuation beam on the left would 


that used in this instance. 
oe absorbed into the fanned portion by two beam coupling, elimunating the split 


appearance of the continued beam. 





*Feinberg shows however, that the electric field may be proportional to the 


gradient of the incident beam intensity [Ref. 27]. 
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To arrive at the configuration of Figure 3.9 a value of An = 0.03 was used. 
MA ieee - 


are characteristic of the value necessary to account for the observed beam fanning in 


max 
= 0.02. Regardless of the profile used, these values of Anna 


BaTiO,. For successful SPPC with the configuration chosen, An pay must be greater 


X 
than 107. Reducing the beam size will obviouslv reduce the necessary value of 


^n but values on the order of 107? to 107? are typical. 


Thus, the simulation has predicted the necessary change in index of refraction 
across the beam path, and the necessary profile of that change. In Chaper V the 
results of the necessary value for An max Will be compared to those possible for the 
crystal. Also, the physical implications of the necessary profile of the index of 
refraction within the beam path, as determined by the simulation, will be discussed. 
The comparison of the information derived from the simulation to that experimentally 
derived from the crystal, will show that the necessary conditions cannot be met in 


BaTiO}. 
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Figure 3.9 PSIM output using the index profile of Figure 3.1. 
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Figure 3.11 Necessary charge distribution and electric field 
to produce the index protile of Figure 3.7. 
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IV. EXPERIMENTS 


A. INTRODUCTION 

In order to investigate beam fanning in BaTiO,, an experiment was designed to 
measure the influence of a change in the incident beam diameter on the fanning 
process. If fanning increases with a decrease in beam diameter, it could be concluded 
that beam fanning is a direct result of the change in index of refraction within the 
beam path, as predicted in Equation 2.25. If, however, the fanning decreases with a 
decrease in incident beam diameter, it could be concluded that beam fanning is actually 
forward-scattering due (probably) to SPS, since, as shown in Equation 2.28, the gain 


per unit length will decrease with a decrease in spot size. 


B. EXPERIMENTAL APPARATUS 

The apparatus used for this work is shown in Figure 4.1. The laser was a Spectra 
Physics Ar’ variable power (20 mW maximum), multimode laser, model #162A-07. 
During the experimentation the laser was operated at all of the visible lines and was 
not optically isolated from the experiment. Upon leaving the laser, the beam was 
directed into a polarization rotator in order to rotate the beam from the ordinary to 
the extraordinary polarization in the crystal. To ensure that only extraordinarily 
polarized light was used, the beam was then directed through a Glan-Thompson 
polarizer. Following the polarizer was a beam splitter, to enable the observation of 
any phase-conjugate reflection, and a variable aperture of the type common in 
photographic lenses. This aperture acted as a spatial filter, reducing the amount of 
scattered light reaching the crystal. The beam was then focused onto the crystal by 
uncoated lenses of either 34 cm or 16 cm focal length. 

The crystal was one of four single crystals of BaTiO, manufactured by Sanders 
Associates and poled into a single domain. Two of the crystals (on loan from the Los 
Alamos National Laboratories) measured approximately 5 x 5 x 5 mm and were 
transparent with a slight yellow tint. Two of the crystals (belonging to the Naval 
Postgraduate School) measured approximately 5 x 5 x 2.5 mm and had no noticable 
ünt to them. The crystal was mounted on a rotatable platform capable of being 
translated 12 cm along the axis of the optical bench on which all of the components 


were mounted. Following the crystal, and connected to the mounting table, was a 
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Knife edge. The knife edge was coated with black tape to ensure that a resonating 
cavity was not formed with the crystal. Behind and to the side of the knife edge, also 
connected to the translation table, was a silicon photodiode connected to an 


electrometer measuring the short circuit current. 


C. EXPERIMENTAL PROCEDURE 

To measure the amount of beam fanning quantitatively, the knife edge was 
aligned in such a way that the beam transmitted through the crystal was blocked from 
the detector. To ensure that the blockage of the transmitted beam was sufficient, the 
table was translated rapidly along its entire range of motion and the detector output 
was monitored to ensure that there was little or no change in the recorded intensity. 
The vibration inherent in the movement and the changing beam diameter of the 
incident beam ensured the lack of significant beam fanning, so that a lack of change in 
the recorded intensity indicated that the transmitted beam was completely blocked 
from the detector for all observed beam diameters. 

Once the crystal, knife edge and detector were in position, the beam was blocked 
and the crystal was flooded with light from a 75 watt incandescent lamp located 
approximately 12 cm away. After approximately one minute the light was turned off 
and the laser beam was unblocked. The intensity of the light reaching the detector was 
then monitored and the peak and steady-state intensities were recorded. Once the 
intensity had reached a steady-state value the beam was again blocked, the crystal was 
translated and flooded with white light, and the procedure was repeated. 

The flooding of the crystal with white light effectively homogenized the charge 
distribution within the crystal and ensured that one reading was not affected by the 
previous one. Special care was taken to ensure that the beam was parallel to the 
crystal table and that all components were mounted on the optical axis of the optical 
bench. The precision translation mechanism on which the crystal table was mounted 
ensured that the beam was incident on the crystal with the same center for all 
measurements. These precautions effectively compensated for any surface 
inconsistencies on the crvstal face except at very small laser beam diameters. 

Both peak and steady-state values of the fanned beam intensity were recorded; 
however, only the peak values were used for final analysis. This decision was made 
because, although the steady state values generally correlated with the peak values, the 
steady state values could be affected by numerous variables. One obvious problem 


with steady-state readings is the possibility of phase conjugation, thus reducing the 
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amount of fanning recorded in the steady state. Also, due to the orientation of the 

detector, 1t was possible (even likely) that the beam would fan across the detector, 

eventually coming to rest past the area it can record; sometimes the fanning was 

observed to exit the face perpendicular to the face the detector was monitoring. To 

reduce all of these effects, the incident beam was originally allowed to enter at only 

small angles to the c-axis (~2° - 5?) and, although alignment was tedious, consistent | 
results were obtained for various crystals, lens focal lengths, laser lines and incident 

angles. Eventually larger incident angles were used (720? - 409?) with similarly 

consistent results. 

To ensure that the observed effect is not influenced by scattering initiated by 
small (e.g. molecular size) crystal inhomogeneities, one configuration was chosen in 
which there was moderate fanning and the fanned intensity was recorded for each of 
the available visible lines of the laser. The result was a mean ratio of incident intensity 
to fanned intensity that was constant across all available wavelengths, thus indicating 
no resonances and no dependence of scattering on wavelength at the available laser 
lines. 

After the data were collected, the beam diameter at each point was measured 
with an Aeronca Electronics Laser Blade and the direct laser intensity was recorded. 
Here, the beam diameter is defined as twice the distance from the center of the beam to 
the point of 10% intensity. For analysis purposes the beam diameter could be divided 
by the cosine of the angle of incidence so that the beam diameter would be that 
actually along the y-axis of the crystal. 

Measurements were also made by placing the detector above the crystal as the 
incident beam diameter was changed. This gave an indication of the total amount of 
scattering, not just asymmetric scattering. In this case, the front surface of the crystal 
was shielded from the detector with black tape to keep specular reflection off of surface 
irregularities from influencing the measurements. Transmission measurements were 
also made by replacing the Knife edge with an aperture, thus allowing only the 


transmitted portion of the incident beam to reach the detector. 


D.. RESULT 

Typical results of the fanning measurements are shown in Figure 4.2. In this 
case the crystal was at an angle of « 5? so the beam diameter need not be adjusted for 
cos§. All intensity recordings have been normalized for comparison purposes. The 


shape of the curves is insensitive to the laser line used, beam angle, choice of crystal, 
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and the divergence or convergence of the incident beam (i.e. whether the crystal is 
placed in front of or behind the focal point of the focusing lens). 

Figure 4.3 is a plot of the intensity seen from above the crystal for various beam 
diameters for converging light. There are two important things to note here. First, 
there is less scattering for smaller beam diameters and second, the sharp discontinuity 
at about three micron beam diameter. This discontinuity is noticable as the crystal 
goes through the focal point of the focusing lens. Figure 4.4 1s a plot of the scattering 
as viewed from above the crystal for diverging light. Note the initial decrease in the 
scattering within the first four microns. This is only apparent with diverging light and 
is seen in all crystals I used. Note also that there is a much broader dynamic range in 
the intensity for diverging light than for converging light. 

Although all of the scattering data are generally insensitive to which crystal is 
used (i.e. the general shape of all the curves for a given situation is the same for all 
crystals), the same cannot be said for the transmitted beam. Figures 4.5 and 4.6 show 
the transmitted intensity for all four crystals plotted together. Crystal Z2 is the only 
one which does not exhibit SPPC and is the only crystal in which the transmitted beam 
intensity increases with beam diameter as a general trend. However, note that all of 
the curves exhibit some inflection point around six micron diameter for diverging light. 
This inflection is not a consistent maximum or minimum; it depends on the crystal. 
Presently I cannot explain the cause of these dissimilar curves. Generally it may be 
noted that the transmitted intensities map the same characteristics as the scattered 
intensity, but inverted, as would be expected. All of these various features will be 
addressed later, but the important aspect of all of the data is obvious; increased beam 
diameter results in increased beam fanning. 

Although these data alone are a strong indication of the lack of index change 
within the beam path due to the photorefractive effect, there exists even stronger 
evidence. By removing the knife edge and detector from behind the crystal, the image 
of the transmitted spot could be observed on a screen approximately three meters 
away. Figure 4.7 shows the image of the spot. which was projected directly onto 
photographic film at different times after the beam was allowed to enter the crystal, 
and prior to the onset of phase conjugation. Figure 4.7a is an image of the beam 
immediately after turn on; notice the spot is generally undistorted. Figure 4.7b is a 
photograph between turn on and steady-state, and Figure 4.7c is a photograph at 
Steady-state. By viewing these images it is clear that the energy of the beam is first 


extracted from the outside of the beam and not from the inside. 
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In addition to the above measurements, the index of refraction of the crystal was 
estimated by measuring Brewester's angle. AtA = 514 nm the published value for the 
extraordinary index of refraction 1s n S 2.424 [Ref. 22], however, my measurements 
consistently resulted in n, = 2.26 0.02. This large discrepancy is still unexplained. It 
is aiso interesting to note that, although there was a distinct minimum at Brewester's 
angle, the reflection did not disappear entirely. The minimum was so slight that it was 
quite difficult to find, and eventually a photodetector had to be used to accurately 
locate it. It 1s possible that both of these observations may be explained by the 
presence of impurities on the surface of the crystal. Specifically, an oil film from the 
hands of the experimenter may be at fault. Additionally, there is the possibility that 
heating at the surface of the crystal due to the high energy density of the incident beam 
may be responsible. Since measurements by others [Ref. 30] of the index of refraction 
of BaTiO, using Brewester's angle have resulted in indices consistent with the reported 
values, either one or both of these explanations may account for the observed 
anomalies. 
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Figure 4.1 Expcrimental Apparatus. 
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Figure 4.2 Dependence of fanning intensity on incident beam diameter. 
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Figure 4.3 Typical dependence of scattered intensity on beam diameter 
as measured from above for converging light. 
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Figure 44 Typical dependence of scattered intensity on beam diameter 
as measured from above for diverging light. 
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Figure 4.6 Transmitted intensity as a function of beam diameter for converging light. 
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igure 4.7a Photograph of projected spot immediately upon turning on the incident beam. 





Figure 4.76 Photograph of projected spot after the onset of beam fanning. 
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igure 4.7c Photograph of projected spot in the steady state. 


V. ANALYSIS 


A. DISCUSSION OF RESULTS 

The simulation in Chapter III and the experiments described in Chapter IV 
clearly show that asymmetric self-defocusing cannot be due solely to a change in index 
of refraction along the beam path within the crystal. Even with the changes of the 
profile for the change in index of refraction proposed in Chapter III, the necessary 
value of An may 18 t00 large according to the published maximum values for BaTiO. 
Inspection of the necessary values of An max from the simulation shows that values on 
the order of 107? are necessary to account for the observed fanning within the crystal; 
however, previously calculated values place an upper bound of ^ 107*, and a value of 
~ 10°? is probably more reasonable? (Ref. 11]. 

This discrepancy may be accounted for only by assuming that the index gradient 
is limited to a very small portion (in some cases as small as "De the size of the input 
beam) of the center of the beam, or by assuming that the observed effect 1s a result of 
another process. [t 1s impossible to justify the first assumption and still achieve the 
necessary index profile, since charge carriers must be liberated from the entire region 
where the brightness of the incident beam exceeds ambient. Therefore, by Laplace's 
equation, the gradient of the electric field within the beam path must be uniform. This 
leaves only the possibility that another process 1s responsible for beam fanning. 

The results of the experiments of Chapter IV reinforce the conclusions made 
from the simulation portion. Figure 4.2 shows that beam fanning decreases with a 
decrease in beam diameter. From Equation 2.25 it is clear that this behavior cannot be 
explained by a change in the index of refraction within the beam path due to the 
photorefractive effect. If this type of change were responsible for beam fanning, the 
the amount of fanning would have to increase with a decrease in beam diameter. 
Additionallv, as shown in Figure 4.7, beam fanning begins on the periphery and ce 
hght is scattered progressing from the ourside of the beam toward the center. Noting 
the index profile as shown in Figure 2.2, this is impossible if the fanning is due to an 
asymmetric change in index of refraction across the beam path, as was originally 
predicted. 


Feinberg does, however, calculate a possible maximum value on the order of 
10° (Ref. 27]. 
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One may argue that all of the observed qualities of beam fanning may be 
explained by assuming “small-scale” self-focusing of the incident beam within the 
crystal. That is, the incident beam may experience a breakup into spots which are of 
diameters small enough to create index gradients of the magnitude shown necessary by 
the simulation in Chapter III. This type of self-focusing has been known for over two 
decades, and is common in materials whose dielectric constant increases with field 
intensity [Ref. 31]. If small-scale self-focusing occurred in BaTiO,, it could explain 
many of the effects that are observed during SPPC (e.g. the presence of small filaments, 
beam fanning, etc.). However, there are serious problems with a small-scale self- 
focusing theory. First, the breakup of optical beams due to self-focusing has been 
shown to be spatially periodic, with the periodicity being dependent upon intensity 
[Ref. 32]. As is demonstrated in Figure 4.7, there is no indication of spatial periodicity 
in the beam as it is transmitted through a BaTiO, crystal. Also, even in good self- 
focusing mediums at very high beam intensities (~ 10’ W/cm’), it is necessary for the 
beam to have a long path length (on the order of tens of centimeters) within the 
medium before self-focusing is observed [Ref. 33]. Even with a very large change in the 
dielectric constant, it is unlikely that self-focusing could occur within a medium of only 
5 mm in length. Finally, there is no evidence in the literature that the dielectric 
constant of BaTiO, changes significantly with the intensity of the incident light 
(although it does change with an intensity gradient, as discussed in Chapter II). It is 
interesting to note however, that the dielectric constant of BaTiO, does change with 
temperature [Ref. 25], and this may cause self-focusing at very high powers and very 
small beam diameters, given a long enough path within the medium. 

From the evidence presented here it follows that asymmetric self-defocusing in 
BaTiO, cannot be due to an asymmetric change in index of refraction along the path 
of the incident beam due to the photorefractive effect. Attributing this phenomenon to 
some other process however, is not so easy. 


B. SCATTERING 

Ecxsssed un Cuapter Ll. the only other theory that is currently available t 
explain asymmetric self-defocusing in BaTiO,, besides the theory just discussed, is the 
theory involving scattering (SPS) (Ref. 24]. Although I have not conclusively 
demonstrated that beam fanning in BaTiO, is due to SPS, there are a number of 


observations which indicate that this type of scattering maybe the origin of the effect. 
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As may be seen in Equation 2.28, if near-forward scattering is the origin of beam 
fanning, then there should be a decrease in fanning with a decrease in incident beam 
diameter, an observation borne out in our experiments. Also, the diversion of energy 
from the incident beam could progress from the periphery toward the center, since 
noise photons originating from outside the beam must interact with the first portion of 
the beam they encounter (Figure 4.7). 

Another interesting and important discovery is that although the shape of the 
fanned intensity curves is similar regardless of convergence or divergence of the 
incident beam (in front of or behind the focal point of the focusing lens), the absolute 
intensity of the scattered light is always observed to be greater for the diverging case. 
This is not perplexing when considering near-forward scattering. One would expect 
more scattering from diverging light than converging because the size of the beam 
increases as it progresses through the crystal for diverging light, thus increasing the 
amount of scattering. Also it is noted that individual intensity readings are not 
consistent, even if there is no no movement of the apparatus between readings. This is 
expected since the scattering process begins with noise photons, and no two situations 
can be exactly alike, even for similar configurations. This inability to exactly reproduce 
any given measurement suggests that the scattering may begin with random index 
changes within the crystal, caused by the stochastic distribution of charge carriers. 

Finally, if stimulated scattering is responsible for beam fanning, then one would 
not expect continuous fanning, but fanning only along paths from which an incident 
photon interacted with the incident beam, thus creating a set of filaments emanating 
from the incident beam. A close examination of any of the photographs of Appendix 
A shows that this 1s indeed the case. Fanning is a set of discrete small beams all 
emanating from the incident beam. We must now discuss curved appearance of the 
dominant-scattered beam during self-pumped phase conjugation, a phenomenon not 
explained by SPS. 

A very close examination of many photographs of SPPC (e.g. Figure A.5) reveals 
that beams which appear to curve often have a series of discontinuities within them, 
Generally, a filament, or set of filaments, leaves the incident beam along the’ direction 
of the c-axis. After a certain distance (usually about halfway to the crystal edge) these 
filaments enter a region where they become fuzzy, and often two become 
indistinguishable. The filaments exit this fuzzy region at an angle different than their 


entrance angle. This type of behavior is explainable by assuming that within this fuzzy 
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region there is stimulated scattering due to SPS. Forward scattered beams enter into a 
region where scattering occurs due to photons arriving from a different angle. The 
beams which continue on at a new angle are the forward scattered beams from this 
second scattering. It is therefore a series of this type of scattering discontinuities which 
gives the characteristic curved appearance of a filament during SPPC. However, this is 
not always the case. 

There exist configurations where this series of discontinuities does not appear to 
account for the curved filaments (i.e. Figure A.6). In a majority of these cases a close 
eXamination reveals that a large filament originates at the incident beam or a series of 
filaments leave the incident beam very close together. The curved appearance comes 
from the coupling of energy to the filament, or portion of the filament, in the direction 
of the c-axis as the filaments progress through the crystal. Therefore the beam doesn't 
actually curve; it simply appears to do so because it is losing energy from one side and 
transferring it to the other. The process is actually a form of self-focusing due to two 
beam coupling. 

I have observed a very few instances where it appears that neither of these 
processes can fully account for the curved appearance of the fanned filaments. In 
these cases one would expect the effect to be due to true beam fanning from the 
change in index of refraction within the filament path. Although the simulation PSIM 
shows that it is impossible for this phenomenon to fully account for the observed 
fanning effect, it would be surprising not to find some evidence of a change in index of 


refraction on the order of 107? across a beam of diameter less than a few microns. 


C. SELF PUMPED PHASE CONJUGATION 

From the discussion of the previous section we may conclude that a change in 
index of refraction within the incident beam path is not responsible for asymmetric self 
defocusing in BaTiO,. It is probable that SPS is responsible for the effect. Therefore 
one may say that the true origin of SPPC in BaTiO, is scattering, probably SPS. The 
cause of the final phase conjugate reflection is not as clear. 

Phase conjugation by stimulated scattering in BaTiO, has been demonstrated 
[Ref. 23], and recently it has been shown how backscattering may dominate forward 
scattering in SPS [Ref. 24]. Additionally, the investigations reported here indicate that 
scattering is responsible for the phenomenon of beam fanning. Therefore one might 
incorrectly conclude that phase conjugation in BaTiO; is due to stimulated scattering. 


With the onset of backward scattering in the crystal there begins a competition of the 
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scattered waves. As the modes begin to build, a mode showing high gain can usurp 
power from one showing less gain. In this way only the modes showing the highest 
gain eventually survive. It is this competition between configurations having high gain 
which is responsible for the amplitude oscillations often seen in a phase conjugate 
reflections (see [Ref. 34,35] or Appendix A). 

If a filament, arising from near forward scattering, has a configuration such that 
it can reflect off the corner of the crystal (i.e. a CAT corner) to create a situation 
where degenerate four wave mixing (DFWM ) is possible in addition to SPS, then one 
would expect that filament to have a gain advantage. This then would be the filament 
to eventually win the competition. So it actually would be expected that both DFWM 
and SPS should contribute to SPPC. 

A review of the conditions of the only published report of phase conjugation due 
solely to SPS [Ref. 23] shows that special precautions were taken to ensure maximum 
coupling. Specifically, an index matching liquid was used which preserved as much of 
the phase conjugate return as possible. Additionally the crystal was surrounded by a 
glass cuvette (containing the index matching liquid) which increased the amount of 
scattered light available for SPS to begin. Even with these advantages, a phase 
conjugate reflection of only ~10% was obtained. This is quite small compared to the 
30% to 50% usually seen from SPPC in BaTiO, in air. This indicates that DFWM 
from the internal corner reflection is the most significant contributor to the phase 


conjugate return observed in SPPC, as was originally proposed [Ref. 2]. 
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VI. CONCLUSIONS 


Experimental evidence does not support the theory that asymmetric self 
defocusing in BaTiO, is a result of an asymmetric change in the index of refraction 
within the path of the incident beam. Instead, forward stimulated scattering, probably 
Stimulated Photorefractive Scattering (SPS), is found to be a more likely explanation. 

Simulation has shown that for an index change to be responsible for the 
phenomenon of beam fanning, the necessary change in the index of refraction across 
the beam must exceed the theoretical maximum, based on the maximum possible space 
charge field. A change in the index of refraction on the order of 10°? is necessary to 
account for the observed effects, whereas the accepted maximum is on the order of 107^ 
£07. 

Experimental results do not support a theory which requires a change in index of 
refraction across the incident beam as the cause of beam fanning. The magnitude of 
fanning is observed to increase with an increase in beam diameter. Also, observations 
of the manner in which the beam fans show that energy is scattered out of the beam 
path beginning with the periphery and not beginning in the center. Neither of these 
observations is consistent with the results expected if a purely photorefractive effect 
were responsible for asymmetric self defocusing. 

Close observations of the process of beam fanning in BaTiO, indicate that beam 
fanning can be explained by a stimulated scattering theory. The presence of filaments 
rather than continuous fanning, scattering beginning with the periphery of the incident 
beam, and discontinuous filaments during phase conjugation are all explicable when 
beam fanning is attributed to stimulated scattering. 

Future experiments in this area should include probing the “fuzzy” regions of the 
observed filaments from above with a low power laser to search for the presence of 
additional index gratings. Also, the actual change in the index of refraction across the 
incident beam path within the crystal remains to be determined experimentally. 
Additionally, further research into the size, type and spacing of scattering centers 
within the crystal is necessary to fully explain the observed difference in scattering of 
diverging and converging light. 
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APPENDIX A 
INTERNAL BEAM PRODUCTION 


l. INTRODUCTION 

The production of internal beams in BaTiO, has, to my knowledge, never been 
systematically investigated. A search of the available literature indicates that the cause 
for this oversight may be that, in general, millimeter sized, collimated beams have been 
used. In all available photographs of BaTiO, used as a SPPCM, the beam has not 
been focused onto the crystal. If the beam is not focused, then the diameter of the 
beam upon entering the crystal is at best the same as when it leaves the laser. It 
appears that when a beam is not focused onto the crystal, auxiliary beams are 
eliminated or obscured by the strong input beam since its diameter is large compared 
to the dimensions of the crystal (sometimes being more than half the size of the 
crystal). However, if the input beam is focused, the beam diameter is small compared 
to the crystal and the immediate effects of the input beam stay localized within the 
crystal, thus allowing the formation and observation of auxiliary beams. There is also 
an immense increase in the energy density of the beam within the crystal upon 
focusing. 

It was using a focused beam that we first discovered the ability of barium 
titanate to support internal auxiliary beam patterns. Extensive investigation by A. V. 
Nowak and myself at Los Alamos National Laboratories indicated that there is a 
plethora of internal beams which may be generated by the nonlinear effects of barium 
titanate during phase conjugation. It appears that phase conjugation is necessary to 
sustain these beams since no auxiliary beams were observed in the absence of phase 
conjugation. In many cases it seems that the beam is initiated by the reflection of the 
phase-conjugate beam off of the face of the crystal as it exits, but this 1s not always 
thule: 

What follows in this chapter is an attempt to report, categorize. and e wp 
many of the observed internal beam patterns and anomalous effects. [t is important to 
note that often it is difficult or impossible to recreate some of the observed effects. It 
appears that not only position and beam profile are important, but also the short term 
history of the crystal. This means that possibly one effect may not be observable 


unless another has (or has not) immediately preceded it. 
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2. EXPERIMENTAL APPARATUS 

All investigations presented in this chapter were performed at Los Alamos 
National Laboratory on an air floated optical bench in free air that was maintained at 
a temperature of 22+ 1°C. The apparatus used is shown schematically in Figure A.1. 

During all experimentation the light source was a Spectra Physics 2020-03 argon 
ion laser, prism-tuned to 514.5 nm and operated in a single longitudinal mode by 
means of an intercavity etalon. The beam at the laser had a diameter of 0.88 mm at 
the !/ ¿ Of peak intensity. 

After leaving the laser, the beam was directed through a half-wave plate that 
served as a beam attenuator when necessary. Following the half-wave plate was an 
apparatus for isolation of the laser from the phase-conjugate beam. The apparatus 
used was a faraday rotator, manufactured by Optics For Research Inc., consisting of a 
25mm long rod of Hoya FR5 glass surrounded by a set of permanent magnets. The 
isolator was set for optimal rotation of the plane of polarization of the beam at 45° 
and placed between two polarizers with axes oriented at 45° to each other. This 
configuration gave a backward attenuation of approximately 500. Without the isolator 
the phase-conjugate beam entered the laser cavity and caused large fluctuations in 
power output as well as mode hopping. The polarization of the beam upon leaving the 
isolator was such that it would enter the crystal extraordinarily polarized, which is 
required for the formation of a SPPCM. 

After passing through the isolator the beam was expanded by a factor of 6.25 
with an expanding telescope. This allowed the insertion of transparencies into the 
beam when desired, enabling the verification that an observed beam is a phase- 
conjugate beam. Also included in the telescope apparatus was a 0.05 mm pinhole that 
could be removed at will. This pinhole provided the ability to make observations with 
a Gaussian beam (pinhole inserted) or with what was eventually termed a mottled 
Gaussian beam (without the pinhole). The mottled Gaussian was so named because, 
although the beam without the pinhole was generally Gaussian, it had a mottled 
appearance. 

Following the telescope, the beam was directed through a beamsplitter with 67% 
transmission at 514 nm. This beamsplitter not only provided the means for viewing the 
phase-conjugate signal as it returned along the incident path, but also provided a 
reference beam that was used to look for frequency shifts of the phase-conjugate beam. 


To look for a frequency shift the phase-conjugate return was reflected from the 
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Figure A.I Experimental Apparatus. 
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beamsplitter and redirected to interfere with a portion of the incident beam reflected 
from the front side of the beamsplitter. The resulting interference pattern was then 
observed for evidence of movement. 

Upon passage through the beamsplitter, the beam, with power of 7545 mW, 
was either modified by a lens, slit, knife-edge or a combination thereof, or simply 
collimated and directed toward the crystal. As the phase-conjugate beam returned 
along the incident path, the beamsplitter reflected a portion of it and directed it to a 
large area silicon PIN photodiode reverse biased at 22 volts and terminated with a 100 
ohm resistor. The output of the detector was sent to a digital millivolt meter, an 
oscilloscope and a chart recorder. When an image was placed in the beam the detector 
was removed and replaced with a white screen for viewing. When frequency shift 
information was desired, the detector was replaced with a mirror and the two beams 
were directed toward a white screen for observation of the resulting interference 
pattern. 

The crystal was any one of four nearly cubic single crystal of BaTiO, measuring 
5+0.2 mm on a side. The sides of the crystals were polished, but it 1s important to 
note that the sides were not exactly the same size and therefore the sides were not 
exactly orthogonal. The crystals were mounted on a platform that could be raised, 
lowered, rotated and tilted. Observations of internal beams were made visually from 
above and were recorded by film or video tape. Both the still and video cameras used a 
f = 35mm Cannon FL lens mounted on a single track bellows extended to provide a 
full frame view of the crystal. 

Many of the photographs presented here were taken from the video tape directly 
off of a television screen using a 6x7 cm format still camera. The desired frame of the 
video was digitized and placed on a hard disk for recall and photographing, the 
numbers displayed in the photographs are a result of the equipment used to do this 
final photography and are not pertinent to the investigation. When this method was 
used, it was possible to make negative images of the desired photographs, which 
reproduce better in print, therefore some of the photographs presented in this work are 
in the negative. There is no other significance to a negative versus positive print. 

In this and in the following sections any reference to phase-conjugate reflectivity 
is the percent return of the beam actually incident on the crystal. The amount of 
specular reflection off of the face of the crystal (which can vary [Ref. 36]) is not 


subtracted out of the calculation. All reflectances have an experimental tolerance of 
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+3%. In the photographs presented, except where noted otherwise, the incident beam 


enters vertically from the top and 1s visible only upon entering the crystal. 


3. RESULTS AND ANALYSIS OF DATA 
a. Experimental variables 

The number of variables involved in SPPC is immense, however I have tried to 
concentrate on those which have emerged as the most important. The most important 
factors in the production of various beams differ, but they. generally may be considered 
as dependent on beam profile or the geometry of the beam-crystal interface. 

Lenses, knife-edges, slits and pinholes were used to change the profile of the 
incident beam prior to the beam entering the crystal. The effect of these modifications 
of the incident beam is especially obvious when the beam profile has a large gradient, 
such as when a knife-edge is inserted into the beam, since this creates a very clear 
delineation in the lighted and unlighted areas of the crystal and introduces new wave 
vectors into the beam. The effect of the lens to crystal distance on phase-conjugate 
reflectivity has been studied by A. V. Nowak and is still not well understood [Ref. 37]. 

The geometry of the beam-crvstal interface is probably the most important of 
all variables in the system. Changes in angle of less than one degree and lateral shifts 
in the beam of less than 0.5 mm can cause drastic differences in internal beams 
configurations. The following examples are presented to indicate the drastic effects 
observed with these type of changes. 

The configuration shown in Figure A.2 1s a typical experimental configuration 
that produced a phase-conjugate signal from the BaTiO, crystal. The phase conjugate 
reflectivity was 33%. Note the fanned beams at the lower left corner of the crystal. 

Figure A.3 demonstrates the extreme effects made possible when changing the 
beam profile by placing a knife-edge in the incident beam. Unless otherwise indicated, 
the knife-edge always intercepted the beam on the right side and was placed behind any 
lenses. Note the extreme curvature of the internal beam observed when blocking half 
the beam with a knife-edge. This also reduced the phase-conjugate reflection to zero. 
At other angles and intersection points. the introduction of a knife-edge produced a 
phase-conjugate retlection where there was none before. 

A large effect is also seen when relocating the beam entrance position by a 
lateral shift in the beam-crystal intersection. As the beam is moved across the face of 
the crystal, the phase-conjugate reflectivity will change, beginning with no phase 


conjugation at one end and ending with the same result at the other. In between the 
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Figure A.2 Normal configuration for BaTiO, exhibiting SPPC. 


two extremes, the phase-conjugate reflectivity may range as high as 60% or be 
nonexistent depending on other factors. Similarly, a change in the angle of the beam 
entering the crystal will have a drastic effect on the phase-conjugate reflectivity and 
beam pattern. It is generally unpredictable and small changes in angle (~4°) have 
been observed to have no effect, to destroy a pattern and eliminate all phase-conjugate 
reflectivity, or to increase the phase-conjugate reflectivity as much as 500%. 

b. Complex Auxiliary Beam Patterns 

~ I. The Total Internal Reflection (TIR) Ring 

One of the most interesting effects discovered was the TIR ring. Figure 

A.4 shows the first example, to mv knowledge, of a totally internally pumped ring 
oscillation within a BaTiO, crystal. This effect was achieved using a Knife-edge to 
obscure half of a mottled Gaussian beam that was then focused through the f= 30 cm 
lens, 29.5 cm from the face of the crystal. The beam was incident at 40° to the axis of 
the crystal and the knife-edge was placed on the side of the beam nearest the crystal. 
The effect was occasionally obtained while obscuring the other half of the beam, and 
was Once obtained with a vertically aligned slit. A. V. Nowak once obtained the ring 
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Figure A.3 Effect of placing a knife edge in the incident beam. 


with just the focusing lens, indicating that the TIR ring may be dependent more on the 
availability of light scattered at a certain angle than on the incident beam profile. 

In the case shown in Figure A.4 the phase-conjugate signal took about 30 
seconds to appear. The beam within the crystal slowly moved to the corner until the 
configuration of the beam arrived as shown, but without the ring. Then, slowly, the 
ring appeared over a three second interval. As the ring appeared, the phase-conjugate 
reflectivity dropped from 76% to 67%. A drop in phase-conjugate reflectivity was 
observed in all examples of ring formation for those cases where the reflectivity was 
measured. In some cases, several nearlv overlapping rings were observed, as in Figure 
A.4. The ring in Figure A.4 was steady, fading occasionally (on the order of once 
every five minutes) but always returning. [t was not overly sensitive to vibration, 
ithough large vibrations would destroy the ring. Rotation of the crvstal destroyed the 
ring but if the rotation was under 1° then the ring would reappear. Translation along 
the crystal on the order of 0.5 mm in either direction did not disturb the ring; further 
translation would destroy it. For other configurations, the ring was only preserved for 


one direction of crystal motion (see Section 4). 
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Figure A.4 The TIR Ring. 


The direction of circulation within the rings was determined by extracting 
some of the light out of the ring into a second crystal of BaTiO,. After producing the 
ring in one crystal, a drop of index matching liquid (methylene iodide (CH,[,) diluted 
1:1 with 2-propanol, n ~ 1.6) was placed onto the surface of a second crystal and the 
second crystal was moved into contact with the first. Since the loss into the second 
crystal was too great to allow the continuation of the ring, the operation was recorded 
by the video camera and the frame of the instant was of contact was studied at a later 
time. Study of the video revealed that the light within the ring could propagate 
counterclockwise, or both wavs simultaneously depending on the particular ring. 
However, in the majority of cases, only the counterclockwise direction was observed. 
Figure A.5 clearly shows this counterclockwise propagation and is exemplary of the 
typical situation found. 

The ring within the crystal is probably produced by two-beam coupling 
between the main beam and scattered light that finds a closed loop via internal 
reflections around the crystal. This somewhat lossy TIR ring takes power from the 


incident beam, which explains the drop in phase-conjugate reflectivity seen immediately 
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Figure A.5 Determination of the direction of the TIR ring 
by extraction of light with another crystal. 


upon its appearance. Of interest also is the difficulty in finding the TIR ring; in some 
crystals the ring pattern could be found in minutes. In others it would take many 
hours. Nowak discovered that in many cases, however, formation of the ring was 
promoted by a slight forward or backward tilt of the crvstal (out of the plane of the 
table). 

The effect of crystal motion on the TIR ring was also investigated by 
Nowak [Ref. 38] and it is interesting to note that the TIR ring appears to seek a square 
Shape. It was observed that the ring will go from rectangular to square upon 
transiation but not from square to rectangular. Also discovered were configurations in 
which the ring exhibited total shape invariance under translation. 

2. The Diamond 

The first of the diamonds, the sharp diamond, is illustrated in Figure A.6 It 
is considerably easier to obtain than the full ring and has been produced with a variety 
ot input beam profiles and orientations. This is thought to be a simple case of two 


beam coupling between the main beam and an internal reflection. Usually when the 
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Figure A.6 The Diamond configuration. 


sharp diamond appears in the crystal the phase-conjugate signal increased considerably 
(up to 50 times what it is Just prior to its appearance). The sharp diamond probably 
begins as a specular reflection of the phase-conjugate signal, which reflects back to the 
incident beam by internal reflection. Since a second point for DFWM 1s created, one 
would expect the observed increase in the phase-conjugate signal upon the appearance 
of the diamond. 

A fascinating and unexplained variation of the sharp diamond is the diffuse 
diamond. The diffuse diamond, shown in Figure A.7, begins in a configuration similar 
to the diamond shown in Figure A.6, but as the crystal is laterally moved, a portion of 
the incident beam composing the sharp diamond stays to create a diamond in the right 
half of the crystal. The diffuse diamond is very weak compared to the other effects 
discussed. but it ts easily visible, especially during its creation. This etfect is similar to 
the motion-invariant effects which will be discussed later. 

3. Internal Specular Reflection 
Often the specular reflection of the incident beam within the crystal 


behaves in unexpected ways. Figure A.8 shows a case where there is visible specular 
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Figure A.7 The Diffuse Diamond configuration. 


reflection of the input beam within the crystal creating a diamond but there is no 
reflection of the secondary beam off of the rear face of the crystal. Absorption of the 
reflected beam by two-beam coupling with the incident beam within the first few 
micrometers of the crystal edge is probably responsible for the lack of specular 
reflection in this and other similar cases. There is a similar theory which accounts for 
the lack of specular reflection of the phase-conjugate beam off of surfaces such as the 
input face of the crystal. Pepper [Ref. 36] proposes that the lack of a reflection of the 
phase-conjugate beam is due to destructive interference between the incident beam 
reflection, and a beam created by the phase conjugation of the reflection of the original 
phase-conjugate beam. This has been borne out in experiment. However, as Figure 
A.9 demonstrates, at certain angles the specular reflection is sufficiently intense to 
make many additional internal reflections. The boundary conditions necessary for the 
lack of a reflection have vet to be fully determined. 
4. Diagonals 
Figure A.10 illustrates the diagonals. The diagonals faded in and out as the 


phase-conjugate reflectivity oscillated from 12% to 21% with a period on the order of 
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Figure A.8 Example of lack of reflection of the incident beam 
Within the crystal. 


5 seconds. The reflectivity was strongest when the diagonals were brightest. When the 
angle of incidence was increased to 56° a similar configuration appeared, but with the 
phase-conjugate reflectivity oscillating from one to 15% and specular reflection off of 
the back face fading as the diagonals became brighter. At another angle of incidence 
in this configuration ( 62?) it was noted that increasing the power of the input beam 
could eliminate all oscillations of the phase-conjugate reflection while eliminating the 
diagonals. This configuration has also been seen without any measurable phase- 
conjugate reflection. 
c. Temporal Effects and Frequency Shifts 

Up to this point this work has considered only examples of the manv internal 
beam patterns that can be obtained by varying the beam profile and orientation. 
Observations of oscillations in the pnase-conjugate amplitude and of frequency shifts 


between the phase conjugate signal and the incident beam will now be addressed. 
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Figure A.9 Example of multiple reflections within the crystal. 


1. Amplitude Oscillations 

Frequently, amplitude oscillations were observed in the phase-conjugate 
signal of many of the previously discussed examples. These include Figures A.2, A.4, 
A.6, A.8, A.9. The oscillations, probably similar to those reported by Valley and 
Dunning in an external resonating cavity [Ref. 39], were most obvious when a mottled 
beam was focused into the sample, but they also appeared in some configurations with 
all beam profiles. The observed oscillations often appeared to be chaotic, and much 
work has been reported recently concerning the observation of chaotic oscillations in 
BaliO, [Ref. 34,35]. However, some of the observed oscillations were very nearly 
periodic (periodic within 10% of some mean period) with periods ranging from 0.01 to 


5 seconds. 


oT 


igure A.I! shows the amplitude of the ohase-conjugate signa! obtained 
with a focused Gaussian beam incident on the center of the crystal 40° to the c-axis. 
These oscillations appeared immediately upon the initiation of phase conjugation. 
These were maintained as shown for a period of about one minute and then abruptly 


changed to the oscillations shown Figure A.12. Within a few minutes the original 
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Figure A.10 The Diagonals configuration. 


oscillations reappeared. Fast oscillations (~10 Hz) amplitude modulated by slower 
oscillations (~0.5 Hz), and oscillations with faster beats superimposed on them were 
also observed. Nowak has shown that the frequency of these amplitude oscillations 1s 
proportional to the intensity of the incident beam (Ref. 37] and this has been confirmed 
by others [Ref. 30]. 

These amplitude oscillations of the phase-conjugate signal are probablv a 
result of frequency shifts of those auxiliary beams that account for the formation of the 
erating. Frequency shifts of the auxiliary beams were proposed bv Feinberg and 
Bacher (Ref. 40] to explain observed frequency sweeping in the output of a ring 
resonator using BaliO, and in an unisolated Ar^ laser producing phase conjugation 
in BaTiO, [Ref. 41]. Thus it is expected that the phase-conjugate reflectivity would 
diminish as the frequency difference between the incident and reflected beams becomes 
larger, since the four wave mixing signal is reduced for large frequency differences as 
shown by MacDonald and Feinberg [Ref. 42]. The nature of and reasons for this 


frequency shifting within the crystal has yet to be explained. 
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Figure A.II Example of quasi-periodic amplitude oscillations 
in BaTiO,. 


2. Frequency Sweeping 

To determine the difference in frequency of the phase conjugate signal, a 
Michelson interferometer arrangement was used. A portion of the incident beam 
interfered with the phase-conjugate signal, as described above, and the resultant 
interference pattern was observed for signs of motion. Normally there was no 
movement of the interference fringes; however, in the case of some configurations 
involving amplitude oscillations (e.g. Figures A.6 and A.10) there was some evidence 
of a frequency difference between the two beams. 

Usually as the phase-conjugate signal began to rise, the interference fringes 
were steady. However, as the signal began to decline the fringes began to move, and 
the more the signal declined the faster the fringes moved, until eventually “heya 
indistinguishable. This sweeping process repeated itself with each oscillation of the 
phase-conjugate reflectivity. The direction of the movement of the interference fringes 
(the direction of frequency shift from the frequency of the incident beam) changed 


arbitrarily, implying that positive and negative frequency shifts were possible. This 
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Figure A.12 Amplitude oscillations which progressed from 
those of Figure A.11. 


correlation of amplitude and frequency shift of the phase-conjugate beam may have 
been explained by MacDonald and Feinberg when they showed that the four-wave 
mixing signal falls off smoothly as the frequency difference between the reflected and 
incident beam becomes larger than ^1 Hz. Although the mechanism for this 
frequency shifting 1s still under investigation, Lam [Ref. 3] has proposed that the bulk 
photovoltaic effect gives rise to an intensity-dependent frequency shift. This implies 
that the frequency is dependent on the amplitude instead of vice versa. Another 
explanation may be that there is consistent movement of the grating when the reflected 
pumpwave forms a loop (see Figure A.2 for example). In this case each circuit around 
the loop will Doppler-shift the pumpwave and the observed effect mav be expected. 
There was only one observed instance where the frequency difference 
between the two beams was non-zero and constant (no photograph available). The 
phase-conjugate return was frequency-shifted from the input beam, but there was no 
sweeping. That is, there was obvious movement of the interference fringes at a 


constant velocity. This effect resulted from a mottled Gaussian beam incident on the 
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crystal at an angle of 35% to the c-axis. However this is an anomalous effect that has 
vet to be repeated or explained. 

3. Influence of Returning the transmitted Beam 

This effect is similar to that reported by Feinberg [Ref. 27] and was 
discovered while trying to stabilize the amplitude oscillations of the phase-conjugate 
signal discussed above. The effect was seen using an arrangement that produced 
nonperiodic amplitude oscillations of the phase-conjugate reflection ranging from zero 
to 46%. In this configuration a 100% reflectance mirror was placed behind the crystal 
directing the transmitted beam back through the crystal, but not directly along the 
incident path. All oscillations stopped, leaving the phase-conjugate reflection steady at 
10%. 

The explanation for this effect is probably found in the competition 
between beam patterns. In any instance, the observed beam pattern is the one which 
exhibits the maximum gain. When oscillation occurs between two competing beam 
patterns it is apparent that, for some reason, the gain of one 1s falling above and then 
below the possible gain of another. The input of a beam of light from any source with 
brightness on the order of the internal beams will cause a change in the gain of the 
system due to some or all of the possible nonlinear effects. In the case of returning the 
transmitted beam, it appears that the gain of one of the competing patterns is 
sufficiently reduced to eliminate it from the possible patterns. 

d. The Free Standing Auxiliary Beam Pattern 

This section presents an occasionally observed result due to a small crystal 
motion. Here the crystal was translated at right angles to the beam path while the 
internal beam pattern was carefully observed for any changes. 

Note the beam configuration seen in Figure A.13. When the crystal was 
shifted slightly to the left, the auxiliary beam pattern moved with the crystal as shown 
in Figure A.14 instead of changing to adapt to the new crystal orientation. Note, 
however, that some slight distortion of the pattern did occur. In this case the phase- 
conjugate reflectivity was originally 32% with small amplitude oscillations on the order 
of 5%. There was no evidence of frequency shifting of the phase-conjugate beam. 
Upon translation of the crystal, the phase-conjugate reflectivity diminished to 23% 
with oscillations of the same amplitude as before the translation. 

After translation, any rotational motion of the crystal destroys the auxiliary 


beam and a new one appears that is similar to the original pattern (i.e. Figure A.13). 
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Figure A.13 BaTiO, prior to translation. 


Also, any translation in the opposite direction, causing the input beam to move over 
the auxiliary beam, will destroy the auxiliary beam and a new pattern will form similar 
to the original. If the input beam is momentarily interrupted and then restored, the 
configuration appearing prior to the interruption returned. 

It 1s important to note that we did not always observe motion-invariance for 
patterns similar to that of Figure A.13. The original pattern as shown in Figure A.13 
was a commonlv occurring one that could be arrived at with any of the beam profiles. 
The angle of the incident beam was not critical and neither was the choice of entrance 
face provided it was parallel to the c-axis. The pattern is usually stable and shows 


phase-conjugate reflectivity as high as 50% at some angles. 
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Figure A.14 Free standing beams created by translating the 
crystal ~0.5mm from the configuration of Figure A.13. 
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APPENDIX B 
PSIM LISTING 


IPSRM S 
! 
{ A SIMULATION OF INTERNAL BEAM RESPONSE IN BARIUM TITANATE 


IBEAM = SIZE OF INPUT BEAM 

IBEAMF = ACTUAL NUMBER OF POINTS IN WAVEFRONT AFTER 
DESIGNATION OF ANGLE PHI 

ISIZE = SIZE OF CRYSTAL 

ICRYS O = CRYSTAL GRID 

KEEP() = ARRAY TO KEEP TRACK OF CURRENTLY OCCUPIED SPACES IN THE 
GRID, 2=COL, 1=ROW 

IPHI = INPUT ANGLE OF INCIDENCE 

PHI = ANGLE OF INCIDENT BEAM (REAL) 

PHIC = ANGLE OF INCIDENT BEAM AFTER ENTERING CRYSTAL 
IX = NUMBER OF COLUMN POINTS ACROSS WAVEFRONT 

IY = NUMBER OF ROW POINTS ACROSS WAVEFRONT 

ADD () = AMOUNT TO ADD TO Y FOR EACH X MOVED 
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OW = ROW NUMBER 
= COLUMN NUMBER 

ITOP = HIGHEST INDEX OF REFRACTION 
IBOT = LOWEST INDEX OF REFRACTION 
GRAD = GRADIENT OF INDEX OF REFRACTION ACROSS THE BEAM 
CRYS = REAL VALUE OF INDEX OF REFRACTION (BEFORE FIXING) 
IMDEX = INDEX OF REFRACTION OF LOCATION 
REALA = PLACE FOR KEEPING TRACK OF AMOUNT TO ADD TO COLUMN 


SOR ROW FOR REALY) WHEN MOVING DOWN THE BEAM 

HIS ENSURES THAT ALL FRACTIONS OF A GRID SPACE 

ARE NOT LOST 

REALY = SEE REALX 

REALAD () = SAME AS REALX ABOVE BUT ONE SPACE FOR EACH POINT. 

USED DURING ACTUAL WAVEFRONT MOVEMENT 

IPERMC() = STARTING COLUMN FOR WAVEFRONT POINT 

IMARKB =MARKER TO KNOW WHEN REFLECTION OFF OF BOTTOM HAS OCCURRED 

IMARKS MARKER TO KNOW WHEN REFLECTION OFF S SIDE HAS OCCURKED 

ICHECK INDEX OF TIME PERIODS 

RSTOP = NUMBER OF TIME PERIODS IN SIMULATION 

IDIR = NUMBER OF TIME PERIODS BETWEEN COMPUTATIONS 

OF DIRECTION OF WAVEFRONT 

MDIR = COUNTER FOR CHECKING AGAINST IDIR 

AX = HORIZONTAL DISTANCE BETWEEN END POINTS OF WAVEFRONT 

YY = VERTICAL DISTANCE BETWEEN END POINTS OF WAVEFRONT 

O SARTRE SEDO TO CHANGE DIRECTION AFTER REFLECTION 
TRAC() = MATRIX USED TO MARK THE INCIDENT BEAM PATH SO THAT 

Dien Lep Seo l UbebyY THE INCIDENT BEAM DOMINATE 
MEAN = MEAN VALUE OF THE INDEX OF REFRACTION 
REF = DX DOR REFRACTION OF THE PHOTOREFRACTOR BEING 


IPLACE = POINT THAT BEAM STARTS IN CRYSTAL 
IGRAPH = MARKER FOR POINTS TO BE PLOTTED (SET FOR 100 POINTS) 
IGRMAX = NUMBER OF ITERATIONS BETWEEN PLOTTING WAVEFRONTS 


CAXIS = FRACTION OF TOTAL INDEX CHANGE DUE TO ANGLE BEAM 
NAKES WET THEE C-AXIS IN THE CRYSTAL 
CINDEX = REAL VALUE OF INDEX OF REFRACTION USED TO COMPUTE 
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THE CHANGE DUE TO THE ANGLE TO C-AXIS 

IOUT() = MATRIX FOR HOLDING VALUES OF PLOTTED POINTS 
ICOUT = COUNTER FOR PLACE OF POINT IN OUTPUT FICE 
ISKIP = MARKER TO CHECK IF WAVEFRONT IS UNDERGOING REFLECTION 
REFLCHKB =MARKER TO KNOW WHEN REFLECTION OF WAVEFRONT STARTS/ENDS 
REFLCHKS = SAME BUT FORSSIUDE 
OP1 = NUMBER OF POINTS IN WAVEFRONT TO LOOK AT FOR RECALCULATION 

THIS EQUALS IBEAMF/PIECE 
ISETMARK = MARKER TO TELL WHEN TO STOP LOOKING FOR REFLECTION 
PIECE = NUMBER OF PIECES OF WAVEFRONT 
ADDY() = ON OF ONE BOX THAT EACH POINT STARTS OUT AHEAD OF 


e-m Com 10 - 4D de GO ÓÉ— a Geom $ bom ¿RD 6-8 Com de es 


INTEGER Ibeam, Ibeamf,Isize,Icrys(500,500) 

INTEGER ie 50 iphi Ix,1y Iro A DN 
INTEGER Ref chkb ( 0), Ikfac,This,That | | 
INTEGER Imdex,Im cu 50),Ipermc(50),Imarkb(50),Itime(50),Iskipb 
INTEGER Imarks(50), ir Moir,Ichnd leout TRI E 1,Piece 
INTEGER Itrac(500,500) ,Mean, Iplace,Igraph, Iout(1,50,101) , Isetmark 
REAL Phi,Phic,Add(50),Grad,Crys,Realx Realad(50 ,Rstop,Grmax 

REAL Xx,Yy,Ref,Caxis(50),Cindex,Rcheck,Kfac,Rea y,Addy(50) 

1 


! 
! NOTE THAT IF BEAM SIZE EXCEEDS 50 YOU MUST REDIMENSION 
! 


ON ERROR GOSUB 5430 ! WRITES FILE TO DISK PITIERROR e 
PRINTER IS 1 II SEITE RIN AUR 
GRAPHICS OFF 
OUTPUT KBD; K; ICLEARST SCREEN 
ON KEY 2 LABEL SAVE &STOP GOSUB 5430 
PRINT TO STOP PROGRAM AND SAVE DATA PRESS f2 


INPUT FACTORS 


6u Gam Dad A A 


INPUT NAME OF DATA FILE TO BE USED ,Name$ 

PRINT DATA FILE: ,Names 

INPUT SIZE OF CRYSTAL TO BE USEDE TISTE 

PRINT CRYSTAL SIZE IS ,ISize ON A STOE 

INPUT BEABMESTZRBOTS _IDeam 

PRINT WITH A BEAM OF SIZE , Ibeam 

INPUT PLACE ON CRYSTAL FOR BEAM TO ENTER IS , Iplace 
PRINT ENTERING AT , Iplace 

INPUT HIGH VALUE OF INDEX IS , Itop 

PRINT LOW VALUE OF INDEX IS ,Ibot : 

PRINT INDICES RANGE FROM pel TO IDOL , 
INPUT INITIAL ANGLE IN DEGREES IS (89 DEG MAX) , Iphi 


TL=TIMEDATE 


! 
! INI TAL IZE VARIABLI 9S 
Rcheck = 0 
Mdir=0 
Ref=2.4 
Icout=0 
Isetmark=0 
Piece=5 ! REM TO CHANGE KK IF YOU CHANGE PIECE 
l 


* 
! 
. 
! 
o 
! 
e. 
! 
o 
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PRINT NUMBER OF PIECES IN WAVEFRONT , Piece 
! 


! 
! 
! PUT FACTORS IN OUTPUT FILE 
! 


Iout(0,0,0)=Isize 
Iout(0,1,0)=Ibeam 
Tout(1,0,0)=Iplace 
Tout(1,1,0)=I Sp 
Iout(0,0,101)=Ibot 
Tout(0,1,101)=Iph1 
Iout(1,0,101)=Plece 


Igraph=0 

Rstop=Itop*2 !SET TIME TOTAL PERIODS 

oa *Isize !USE TWO LINES TO KEEP FROM INTEGER PROBLEMS 
PRINT TOTAL TIME = , Rstop 

Grmax=Rstop/100 ! SET FOR 100 WAVEFRONTS PLOTTED 

GOSUB "5170" ISET UP GRAPHICS DISPLAY 

! 


! TIME BETWEEN NONLINEAR CALCULATIONS 
eee cance? 


! 
INITIALIZE ARRAY 


FOR I=1 TO Ibeam 


EE PHI FROM DEGREES TO RADIANS AND FIND ANGLE INSIDE CRYSTAL 
i=Iphi 

Phi=Phix (24PI 360) 
Phic=(SIN(Phi)/Ref) 
Phic=ASN(Phic 


! 

! CHANGE SIZE OF BEAM DEPENDENT ON THE INPUT ANGLE 

4 

e e nic 

Y=Ibeam*SIN(Phic 

Iy=Y 

Ix=X 

Ibeamí=Ix_ 

SS !IISIZE OF PIECE TO LOOK AT WHEN CALCULATING DIRECTION 
Addi-Add(1) 

TE Iplace>(Isize-Ibeamf) THEN ! MAKE SURE IPLACE IS NOT TOO LARGE 

ies =Isize-Ibeamf 
END IF 
1 


: 
! 
! 
! 

Irow=0 
Realy=0 
Icol-Iplace 


FILE WAVEFRONT COORD 
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FOR I=1 TO Ibeamf 
Keep(2,1I)-Icol 
Rea c +Add(1) 
Add(1)=Add(1) 
IF Realy>Irow+1 THEN Irow=Irowtl 
A eo eee 
Keep(1,I)=Irow 
Icol=Icol+1 
B DI 


SET UP CRYSTAL TO MEAN INDEX 


ean=ABS ((Itop+Ibot)/2) 
FOR I-1 TO Isize 
FOR J=1 TO Isize 
Icrys AN 
Itrac(1,J)=0 


N 
! 
i 
] 
I 
| AAKAKSET UP INDEX OF REFRACTION FOR PATH OF INCIDENT BEAM **** 
! 
! 
! 
! 
M 


! 
Grad-(Itop-Ibot)/2 !INDEX OF REFRACTION GRADIENT 
FOR I-1 TO IDeamf 
Irow=Keep(1,I 
Icol-Keep(2,I 
Crys-Grad^SIN(((Ibeamf£/2)-(I-1))*(PI/Ibeamf)) 
Crys=Crys+Mean 
i A ery I DIE 


SET UP VALUES DOWN THE BEAM PATH 


OR I=1 TO Ibeamf 
pm E 
Irow-Keep(1,I 
Imdex=Icrys(Irow,Icol) 
Imdexp (I )=Imdex 
Itemp=Icol 
Realx=0 
FOR J=Irow TO Isize 

leuc p Item. -Imdex 
Itrac I r =1 
Realx=Realx+Add(I) 
T , 
IF Itemp»Isize THEN Itemp-Isize 
IF Itemp«O THEN ItempzO 
NBI ere) 
SUEDE 


ps — sa one owes ones 2% 
[rj 


SET THE STARTING COLUMNS FOR THE WAVEFRONT AND 
INITIALIZE FLAGS 


FX] ome o-a o-a ema o-a Ome ome 


OR I=1 TO IBEAMF 
ipeum EE 
Imarkb(1)=0 
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Imarks(1)=0 
Tus 


XX*XX*kBEGIN MOVEMENT XXARXARARARAAAAAARA RARA KA AAA KK AR 


MAKE SURE BEAM IS STILL IN CRYSTAL 


FOR I=1 TO Ibeamf 
Ir Imarkb(1)=1 
IF Keep(l, 1)<2 THEN GOTO 5140 
END IF 
IF Imarks(I)=1 THEN 
IF Keep(2,1I)>Isize-10 THEN GOTO 5140 
END IF 
ae T 


1 AARKÁXÁMOVEMENT XAARARARARARA RARA RARA RARARR ARA RARA AA RARA A 


l 
! 


WHILE Rcheck«Rstop 

! 

FOR I=1 TO Ibeamf s 
Icol= E I 


Irow=Keep(l 
pOr y (bos Icol) 


MARK REFLECTION IF NECESSARY 
IF Icol<=1 THEN 


IF Imarks(1)=0 THEN 
Imarks 


$e t-a den ¿eo 


NONESBES TT IS TIME 


IF Itime(1)<Imdex THEN 
E =Itime(I)+1 


Itime(I)=0 

Realad(I) pos cade) 

IF Imarks(1)=0 THEN 
Icol=I ier peated) 
PE zoe «0 THEN Icol=0 


ELSE 
Icol=Realad(I) 

END IF 

Keep(2,10=Icol 

IF markab(1)- O THEN 
E (1,1)=Irow+1 


E 
Keep(1,I)=Irow-1 
IF 


-p Gam dr m 
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PUT POINTS IN FILE AND DISPLAY GRAPHICS 
EE Aris TINE TOTORO 


&— er Qo e SuM 


FS raph+1 
D raph Grmax THEN 
E cout+l 
FOR L=1 TO Ibeamf 


Igraph=0 
Pede Gee LeO ZE me !X COORD OF POINT 
al Icout)=Isize-Keep(1,L) !Y COORD OF POINT 
PLOT Iout(0',L,Icout), Iout(1 ,L, Icom 

NEXT L 


EU TE 


! 

l 

! 

! 

e 

! CHECK IF TIME TO RECOMPUTE DIRECTION **XXXXXXXXAXAXXXX XX 


MDIR=Mdir+1. 
IF depa THEN 


Mdir=0 
PRINT TABXY(1,11), CURRENT TIME PERIOD = , Recheck 


OR K=0 TO Piece-1 
kkkkkkkkkkkkkkkk SET FOR 5 PIECES žžkkkkkkkkkkkkkkk 


ome ome ones PY] ome cme om 0-0. 0. 0-. 


IF K=0 THEN Kk=0 
IF K=1 THEN Kk=4 
IF K=2 THEN Kk=1 
IF K=3 THEN Kk=3 
IF K=4 THEN Kk=2 


!IF K=9 THEN Kk= 
This-Kk*o e 
That=(Kk+1)*0 
IF Kk=Piece-1 THEN That=Ibeamf 


FIRST CHECK TO MAKE SURE THE FRONT IS NOT PARTIALLY REFLECTED 


- omm 06-a 0 1 


Iskipb=0 
Isetmark=0 
IF Reflchkb(K)20 THEN 
FOR M=This TO That 
IF Imarkb(M)=1 THEN 
Isetmark=1 
FOR N=This TO That 
IF Imarkb(N)=0 THEN Iskipb=1 
NEXT N 
END IF 
NEXT M 
END IF 
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IF ERE AND Iskipb=0 THEN Reflchkb(K)=1 
Iskips=0 

Isetmark=0 

IF Reflchks(K)=0 THEN 

FOR M=This TO That 

IF Imarks(M)=1 THEN 


FOR N=This TO That 


Isetmark=1 


IF Imarks(N)=0 THEN Iskips=1 
EXT N 


END IF 
NEXT M 


mai y Gow 6 


! 


IF Isetmark-1 AND Iskips=0 THEN Reflchks(K)=1!SKIP THIS IF NOT NECESSARY 
! 
END IF 


e dl e de 


e—a = a Out Omat de 


IF Iskips=0 AND Iskipb=0 THEN !IF WAVEFRONT IS NOT PARTIALLY REFLECTED 
SLOPE OF PIECE OF WAVEFRONT IS YY/XX 


viz =Itime Bee eee Een ST SUE 
Y2=Itime(That 1,That) ,Keep(2,That 
Yy=Yy+Y1-Y2 

IF Imarks(This)=0 

"t E T G "Realad(That))- (Ipermc(This)-Realad(This)) 

Xx-Realad(That)-Realad(This) 

END IF 

IF Xx>1 THEN 

FOR G=This TO T 
AAA (C) ABS (Yy/ 2x) !IFIGURE NEW DIRECTION 
Caxis(G)-7Xx/Opl 

IF Kk=Piece-1 THEN Caxis(G)=Xx/(Ibeamf£-(K*0p1)) 


NEXT G 
PRINT Add(G-1) 
END IF 


reine is)/torys{Keep(t 3 (Keep(1, E 
í 


/Icrys (Keep 


| 
! 
] 
! 
FOR I-This TO That 
sl H 
Irow-Keep(l,I 
Imdex=Imdexp (1) 


! 
Ltemp=Icol 
Realx=0 


CHECK FOR REFLECTION 


IF Imarkb(I)=0 THEN 
FOR J=Irow TO Isize 
Icrys(J,Itemp)=Imdex 
Rea x=Realx+Add(I) 
PESE =Icol-Realx 
IF Itemp<0 THEN Itemp=0 
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IF Itemp>Isize THEN Itemp=Isize 
Itrac(J,Itemp)=1 


NEAT J 
! 
! 
ELSE ! INCIDENT BEAM DOMINATES 
FOR J=0 TO Irow 
Ichng-Irow-J 
IF dc erence )=0 THEN 
Icrys Ichng Pte ES 
Icrys(Ichng,Itemp)=Icrys(Ichng,Itemp) 
ma EF 
Realx-Realx-Add(I) 
Itemp Ico ira 
IF Itemp<0 THEN P 
IF Itemp»Isize THEN Itemp-Isize 
END IF 
NEXT I 
END IF 
NEXT K 


! CHECK FOR END OF RUN XXXXXXXXAXXXHAXAXAXA AX 
END WHILE 
T2-TIMEDATE 

PRINT TOTAL RUN TIME = ,T2-Tl 
GOSUE 5360 Waite, FILE TO DISK 
T 

! SUBROUTINE TO INITIALIZE GRAPHICS 

GRAPHICS ON 

GCLEAR 

SHOW 0,Isize,0,Isize 

PLOT 0) O02 

FOR Q-0 TO Isize 

PLOT Q,0,*] 

NEXT 

FOR Z-0 TO Isize 

PLOT Isize,Z,-1 

NEXT Z 

FOR Q=Isize TO O STEP -1 

PLOT Q, Isize,-1 

NEXT Q 

FOR Z=Isize TO O STEP -1 

PLOT D ou 


! SUBROUTINE TO OUTPUT FILE TO DISK 

MASS STORAGE IS :,700,0 

CREATE BDAT Name$,10404,2 

ASSIGN Path_1 TO Names 

OULEUD Sr ah our Gs} 

ASST eaten ie LO ^ 

RETURN 

! SUBROUTINE TO OUTPUT FILE IF ERROR OCCURS 
MASS STORAGE [5 70070 

CREATE BDAT Names ,10404,2 

ASSIGN Path_1 TO Names 

OUTPUT Path_1;Iout(*) 

ASSIGN Path_1 TO * 

PRINT ERROR! PROGRAM TERMINATED 

T2=TIMEDATE 

PRINT TOTAL RUN TIME = ,T2-T1 . 
STOP 
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5520 
5530 


RETURN 
end 
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